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Abstract 

The recently developed on-shell bootstrap for computing one-loop amplitudes in non- 
supersymmetric theories such as QCD combines the unitarity method with loop-level on-shell 
recursion. For generic helicity configurations, the recursion relations may involve undetermined 
contributions from non-standard complex singularities or from large values of the shift parameter. 
Here we develop a strategy for sidestepping difficulties through the use of pairs of recursion rela- 
tions. To illustrate the strategy, we present sets of recursion relations needed for obtaining n-gluon 
amplitudes in QCD. We give a recursive solution for the one-loop n-gluon QCD amplitudes with 
three or four color- adjacent gluons of negative helicity and the remaining ones of positive helicity. 
We provide an explicit analytic formula for the QCD amplitude Aq-i{1~ , 2~, 3~, 4^, 5^, 6^), as well 
as numerical results for A7;i(l-, 2-,3-,4+, 5+, 6+, 7+), yl8;i(l", 2-,3-,4+, 5+, 6+, 7+, 8+), and 
As-i{l^ , 2^, 3^, 4~, 5^, 6^, 7^, 8^). We expect the on-shell bootstrap approach to have widespread 
applications to phenomenological studies at colliders. 

PACS numbers: 11.15.Bt, 11.25.Db, 11.25.Tq, 11.55.Bq, 12.38.Bx 



I. INTRODUCTION 

The success of the forthcoming experimental program at CERN's Large Hadron Colhder 
will depend in part on theoretical tools. Our ability to find and understand new physics at 
the TeV scale will rely on the quality of predictions for a variety of known-physics processes. 
A classic example is the H^ + 4 jet background to top-quark production. Tools to perform 
higher-order corrections to a wide variety of processes in the component gauge theories of the 
Standard Model will play an important role. Tree-level scattering amplitudes provide the 
basic predictions for cross sections for Standard Model processes. However, next-to-leading 
order (NLO) QCD corrections are typically quite large. One-loop QCD amplitudes, which 
enter at NLO, are therefore needed in order to reduce theoretical uncertainties to the level 
of 10% or so. An important set of Standard Model backgrounds to new physics dictates the 
computation of new one-loop amplitudes for processes containing one or more vector bosons 
(Ws, Zs, and photons) and multiple jets. 

Experience has shown that while methods relying on direct analytical evaluation of Feyn- 
man diagrams can be used for five-point processes, they have not proven powerful enough to 
compute six-point processes or beyond in QCD. The recent development of semi-numerical 
approaches [1, 2, 3] shows promise for improving traditional capabilities. All helicity con- 
figurations for the six-gluon amplitude have been evaluated numerically in this way, and 
numerical results presented for a single phase-space point [3]. These results are also of util- 
ity in confirming analytic expressions. (For other numerical or semi-numerical approaches, 
see ref. [4].) 

On-shell methods for computing amplitudes can be much more efficient than Feynman 
diagrams, because they avoid gauge non-invariant intermediate states and instead focus 
on the key analytic properties that any physical amplitude must satisfy. The unitarity- 
based method [5, 6, 7, 8] was applied long ago, not only to six-point processes, but also 
to all-multiplicity amplitudes, for particular configurations of external helicities. Early ap- 
plications of the method were generally restricted, for practical reasons, to supersymmetric 
theories or to the polylogarithmic part of QCD amplitudes. This practical restriction arose 
from the greater complexity of D-dimensional unitarity calculations, required for full QCD 
amplitudes in this approach. 

A key feature of the unitarity method is that new amplitudes are constructed with only on- 



shell tree-level amplitudes (which are generally quite simple) as inputs. A number of related 
techniques have emerged in the past two years, including the application of maximally- 
helicity-violating (MHV) vertices [9, 10] to loop calculations [11, 12] and the use [13, 14] of 
the holomorphic anomaly [15] to evaluate the cuts. 

More recent improvements to the unitarity method [16, 17, 18, 19] use complex momenta 
within generalized unitarity [20, 21, 22], allowing, for example, a simple and purely alge- 
braic determination of all box integral coefficients. (The name 'generalized unitarity', as 
applied to amplitudes for massive particles, can be traced back to ref. [23].) In ref. [18], 
Britto, Buchbinder, Cachazo and Feng developed efficient techniques for evaluating generic 
one-loop unitarity cuts, by using spinor variables and performing the cut integration via 
residue extraction. Quite recently, Britto, Feng and Mastrolia [19] further extended these 
techniques and completed the computation of all cut-containing terms for the six-gluon he- 
licity amplitudes. The cut-containing terms for other helicity configurations, and for other 
components of the amplitudes, were obtained in refs. [5, 6, 18, 24, 25, 26]. The only terms 
now missing in the analytic expressions for the six-gluon amplitudes are the pure-rational 
ones. The computation of the rational terms, in these and more general amplitudes, is the 
subject of this paper. 

In a previous paper [27], three of the authors presented a systematic, recursive bootstrap 
approach to making high-multiplicity QCD calculations practical within the framework of 
the unitarity-based method. It complements the use of four-dimensional unitarity for loga- 
rithmic and polylogarithmic terms with an on-shell recursion relation [28, 29, 30, 31] for the 
purely-rational terms. This approach systematizes a unitarity-factorization bootstrap previ- 
ously applied to the amplitudes for e+e~ — > 4 partons [20]. It has already been used to solve 
for infinite sequences of one-loop n-gluon helicity amplitudes, in particular the MHV am- 
plitudes containing two color-adjacent negative-helicity gluons and (n — 2) positive- helicity 
ones [32]. These papers do not explain how to attack more general helicity configurations. 
That is the purpose of the present paper: to extend the range of applicability of the recursive 
bootstrap method to cover as generic a helicity configuration as possible. 

Recursion relations have long been used in QCD [33, 34], and are an elegant and effi- 
cient means for computing tree- level amplitudes. Other related approaches [35], as well as 
computer-driven approaches such as MadGraph [36], have also been employed. Stimulated 
by the compact forms of seven- and higher-point tree amplitudes [22, 37, 38] that emerged 



from studying infrared consistency equations [39] for one-loop amplitudes (computed using 
the unitarity-based method), Britto, Cachazo and Feng wrote down [28] a new set of tree- 
level recursion relations. The new recursion relations differ in that they employ only on-shell 
amplitudes (at complex values of the external momenta). A simple and very general proof 
of the relations, using special continuations (shifts) of the external momenta in terms of a 
complex variable z, was then given by Britto, Cachazo, Feng and Witten [29]. The power of 
this type of recursion relation follows from the generality of the proof, which relies only on 
factorization and Cauchy's theorem. (The numerical efficiency of these recursion relations, 
with respect to the older, off-shell recursion relations [33, 34] and those based on MHV ver- 
tices [9, 40], has been studied recently [41].) On-shell recursive methods have also yielded 
compact expressions for tree amplitudes in gravity [42] as well as gauge theory [43] , and have 
been extended to theories with massive scalars and fermions [44, 45]. They even provide a 
derivation [10] of the Cachazo-Svrcek-Witten representation of amplitudes in terms of MHV 
vertices [9]. Many of these developments, as well as the resurgence of interest in unitarity 
methods, were inspired by the development of twistor string theory [46]. 

The unitarity-based method [5, 6] turns a general property of field theories — the unitarity 
of the (perturbative) S'-matrix — into a practical technique for computing cut-containing 
terms in amplitudes. In a similar spirit, on-shell recursion relations turn another general 
property — factorization on poles in intermediate states — into a technique for computing 
rational terms in amplitudes. The idea of using factorization as a computational tool goes 
back to the computation of the Z ^ A parton one-loop matrix elements [20] (or equivalently, 
by crossing, the virtual diagrams for pp -^ W, Z+jets), wherein all terms consistent with the 
helicity assignments were written down, and collinear limits used to isolate the correct ones 
and their coefficients. This approach gets harder to apply as the number of external legs 
increases, because of the difficulty of finding terms with the correct factorization properties. 
The one-loop on-shell recursion relations [27] provide a practical and systematic method for 
constructing the rational terms, avoiding this difficulty. Moreover, in special cases, when 
certain criteria on the unitarity cuts are satisfied [26] , it is also possible to obtain the rational 
coefficients of the cut-containing (poly) logarithmic terms via on-shell recursion relations. 

The factorization properties of one-loop amplitudes in gauge theories, as a function of 
real Minkowski four-momenta, have been known for a long time. We may distinguish two 
different cases. In the first case, dubbed "multi-particle" factorization, the momentum 



going on shell is a sum of three or more external momenta. In the second case, called 
"collinear" factorization, it is a sum of two momenta. The standard derivations describe 
how amplitudes factorize in either of these limits, when all momenta involved are real. The 
implementation of on-shell recursion relations requires a generalization of these factorizations 
to complex momenta. This generalization is straightforward, both at tree level and at one 
loop, for multi-particle factorization. The generalization is also straightforward at tree level 
for collinear factorization. This is no longer true at one loop. 

The heuristic reason why collinear factorization is more intricate with complex momenta is 
that one cannot define a nonsingular all-massless three-point kinematics with real momenta, 
while one can using complex momenta. At the loop level, the complexity of complex collinear 
factorization is refiected in the appearance of double poles and 'unreal' poles in scattering 
amplitudes [30, 31]. As yet, we have no general theorems providing universal factorization 
formulae in these "non-standard" cases. In previous computations [27, 32] of one-loop ampli- 
tudes with two color-adjacent negative legs, these problems could be sidestepped by making 
special choices in constructing recursion relations. Within the framework of ref. [27], we 
must choose momentum shifts under which the amplitude vanishes at large shift parameter 
z. Otherwise, the contour integral over z which gives rise to the recursion relation would 
receive an undetermined contribution from large z. In the case of MHV n-gluon amplitudes, 
which contain only two negative-helicity gluons, it is possible to make such a choice and yet 
avoid channels with unknown factorizations [27, 32, 47]. 

For general helicity configurations, however, it is no longer possible to do this. It might 
seem that we should therefore study the 'difficult' channels, and attempt to derive a universal 
form for their complex- momentum factorization. It turns out, however, that it is easier to 
relax the other requirement, that of a vanishing amplitude at large shift parameter z. Indeed, 
in ref. [30] it was shown that, if we somehow knew the large-z behavior of an amplitude in 
a recursion, then a non-vanishing behavior posed no problems; the recursion relations still 
reconstructed the remaining terms in an amplitude correctly. Our aim here is to show how 
to determine the large-z behavior of amplitudes from scratch. We will do so by using an 
auxiliary recursion relation, constructed by considering pairs of momentum shifts, one in the 
parameter z and a second involving different external legs and another parameter w. With 
these additional terms in hand, we can follow the approach of ref. [27] for the remainder 
of the calculation, computing recursive and overlap diagrams to add to the cut-containing 



terms. 

As an illustration of our method, we will compute one-loop corrections to a class of 
next-to-maximally-helicity-violating (NMHV) n-gluon amplitudes in QCD, those with three 
adjacent negative helicities in the color ordering, A„.j^°°^( 1-+ ■ ■ ■ ++)• Under a super- 
symmetric decomposition [48], these amplitudes may be thought of as composed of A/" = 4 
and A/" = 1 supersymmetric pieces together with a non-supersymmetric (A/" = 0) scalar loop 
contribution. The A/" = 4 contributions were computed in ref. [37], and the A/" = 1 terms 
in ref. [49]. The logarithmic parts of the A/" = scalar loop amplitudes were determined 
in ref. [26], by constructing an on-shell recursion relation for integral-function coefficients 
appearing in the amplitudes. We shall complete the QCD computation in this paper by 
obtaining the rational- function contributions. 

We also describe a recursive solution for the rational-function parts of the scalar loop am- 
plitudes with four color- adjacent negative helicities, using the logarithmic terms computed 
in ref. [26] as a starting point. We have computed in this way the A/" = terms in the 
eight-gluon amplitude Aq\°°^{ 1-+++). 

We present numerical values for the six-, seven- and eight-gluon amplitudes with "split" 
helicity configurations, in which all the negative helicities are color-adjacent, as a reference 
point for future implementations of these amplitudes in phenomenological studies. 

This paper is organized as follows. In the next section, we review notation and the 
organization of color-ordered amplitudes used in this paper. In section III, we present a 
known five-point amplitude, to illustrate and guide our strategy for obtaining the rational 
parts of one-loop amplitudes with general helicity configurations. In section IV, we then 
apply this strategy to determine a sample six-point amplitude. Before continuing to more 
general cases in section V, we review and extend the on-shell bootstrap formalism [27] to 
cases where the shifted amplitudes An{z) do not vanish for large shift parameter z. In 
section VI, we observe various empirical properties, which we use to construct a procedure 
for general helicities, focusing on n-gluon amplitudes. As a non-trivial confirmation of the 
general procedure, in section VII we present examples of applications of our procedure 
for determining the behavior of amplitudes for large values of the shift parameter. This 
procedure is then used in section VIII to determine a recursive solution of the rational 
functions for n-point amplitudes with three nearest-neighboring negative helicities in the 
color ordering. We also describe a recursive solution to the eight-gluon amplitude with 



four color-adjacent negative helicities, Ag.]^°°^( 1-+++). In section IX we present 

numerical values of the scattering amplitudes at select kinematic points. In section X we 
present our conclusions and outlook for the future. We include an appendix collecting 
previously computed amplitudes that feed into our recursive computations. 

II. NOTATION 

In this section we summarize the notation used in the remainder of the paper. Following 
the notation of previous papers [27, 30, 31], we use the spinor helicity formalism [50, 51], in 
which the amplitudes are expressed in terms of spinor inner-products, 

(j7) = {r\i^) = u4k^)u4k), [ji] = {j+\n = u4k,)u4k), (2.1) 

where u±{k) is a massless Weyl spinor with momentum k and positive or negative chirality. 
The notation used here follows the QCD literature, with [ij] = sign{k^k^) {j i)* so that, 

{tj)[jt] = 2ki-kj = s,j. (2.2) 

Our convention is that all legs are outgoing. We also define, 

Xi = u+(ki), \i = u_{ki). (2.3) 

We denote the sums of cyclicly-consecutive external momenta by 

Kl., ^ K + kl, + --- + k^-, + k^ , (2.4) 

where all indices are mod n for an n-gluon amplitude. The invariant mass of this vector is 

s^.., = Kl, ■ (2.5) 

Special cases include the two- and three-particle invariant masses, which are denoted by 

Sij ^= -tiij ^= (fcj + kj) = 2fcj ■ Kj, Sijk ^= [ki + Kj + K}^) . (2.0) 

We also define spinor strings, 

{r\{a±b)\r) = {ta)[aj]±{tb)[bj] , 
{i+\{a + b){c + d)\j-) = [ia]{a-\{c + d)\j-) + [ib]{b-\{c + d)\r) . (2.7) 



We use the trace-based color decomposition of amplitudes [51, 52, 53, 54]. For tree-level 
amplitudes with n external gluons, this decomposition is, 

AT{{KKa,})=f-^ Y^ Tr(T"-W---TMn)X^<=(a(l^\...,n'^")). (2.8) 

Here g is the QCD coupling, Sn/ Zn is the group of non-cyclic permutations on n symbols, and 
jhj denotes the j**^ gluon, with momentum /cj, helicity hj, and adjoint color index a^. The 
SU(A'^c) color matrices in the fundamental representation are normalized by Tr (T"T^) = 5°"^. 
For spin-J adjoint particles circulating in the loop, the color decomposition for one-loop 
n-gluon amplitudes is given by [55], 

\n/2\+l 

At^°'^\{KKa,})=g-Y.nj J^ Y. Gr,,c(a) 45^ • (2-9) 

J C=l crdSn/ Sn;c 

The notation in eq. (2.9) has been described repeatedly elsewhere [27, 30, 31, 55]. 
Here we just note that we need to compute only the leading-color partial amplitudes 
A^.[(l''% . . . ,n^"), because the subleading-color partial amplitudes for a gluon in the loop, 
An;c for c > 1, are given by a sum over permutations of the leading-color ones [5]. The 
analog of eq. (2.9) for fundamental-representation particles in the loop (such as quarks, with 
spin J = 1/2) is also expressed in terms of Al^.[, 

A'r''{{h,K,a,})=g^ Y ""J Yl Tr(TMi)...TM.))4f^(a). (2.10) 

J=0,l/2 aeSn/Zn 

The contributions of different spin states can be rewritten in terms of supersymmetric 
and non-supersymmetric parts [48], 

<1 =<r-4<f + <r°- (2.12) 

The non-supersymmetric amplitudes, denoted by A/" = 0, are just the contributions of 
a complex scalar circulating in the loop, A^.^^ = Aj^. The supersymmetric and non- 
supersymmetric pieces have different analytic properties. The supersymmetric pieces can be 
constructed completely from four- dimensional unitarity cuts [5, 6] and have no additional ra- 
tional contributions. The polylogarithms and logarithms of the A/" = non-supersymmetric 
contributions may also be computed from the four- dimensional unitarity cuts. (In certain 
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cases, the coefficients of integral functions containing the logarithms and polylogarithms 
may instead be determined recursively [26].) 

The leading-color QCD amplitudes are expressible in terms of the different supersym- 
metric components via, 

^„^r = <r' - 4<r^ + (1 - ^^i^Xr" + ^ (<r^ - <r°) + ^<r° , (2.13) 

where Uf is the number of active quark flavors in QCD. We also allow for a term proportional 
to the number of active fundamental representation scalars n^, which vanishes in QCD. We 
regulate the infrared and ultraviolet divergences of one-loop amplitudes dimensionally. (In 
this paper, we will not treat divergent and finite parts separately.) The regularization- 
scheme-dependent parameter Sr specifies the number of helicity states of internal gluons to 
be (2 — eSpi). For the 't Hooft-Veltman scheme [56] 6^ = 1, while in the four- dimensional 
hehcity (FDH) scheme [6, 57, 58] 5r = 0. 

The amplitudes in this paper have not been renormalized. To perform an MS renormal- 
ization, subtract from the leading-color partial amplitudes A^.^ the quantity, 

n — 21 f 11 2 Hf 1 rig 



cr 
with the universal prefactor 



3Nr 3N, 



AT\ (2.14) 



III. A FIVE-POINT AMPLITUDE THE HARD WAY 

A. Overview of the On-shell Bootstrap 

In this paper we will continue the development of the method of ref. [27] for obtaining 
complete one-loop amplitudes in QCD and other non-supersymmetric theories. In refs. [27, 
32], amplitudes with two color- adjacent negative-helicity legs were considered. In that case, 
it was possible to choose a shift so that: 

1. the recursion relations do not contain any terms with non-standard complex factor- 
izations, and 

2. the amplitude vanishes for large values of the complex shift parameter z. 

10 



However, for general helicity configurations it is not possible to satisfy both conditions at 
once. Here we will provide a simple procedure for sidestepping this apparent difficulty. 

Before we attempt to calculate some six- and higher-point helicity amplitudes with three 
negative helicities, for which both conditions cannot be satisfied, let us examine a five-point 
amplitude, A^.^ (1^,2^,3^,4^,5"*"). We shall focus on the scalar-loop contribution, A^j"^, 
which is the only component not computable from the D = 4 unitarity cuts. We know the 
answer for this amplitude [48], which makes it a good test case. Of course, it is maximally- 
helicity- violating in the conjugate spinors, and therefore could be computed as in ref. [27]. 
Here, we will instead compute it in a way that foreshadows our computation of higher-point 
amplitudes. 

We now briefiy review the construction of ref. [27]; in section V we will present a more 
systematic review and extension of the on-shell bootstrap. We first choose a complex-valued 
shift [29] of the momenta of a pair of external particles, kj -^ kj{z), ki -^ ki{z). We describe 
a [j, /) shift in terms of the spinor variables A and A defined in eq. (2.3), 

[j, I) : \j -^ \j - z\i , \i^ \i + z\j . (3.1) 

To compute the amplitude, we must first examine the cut terms, obtained from the 
unitarity method and/or recursion on the coefficients of integral functions, for the presence 
of spurious singularities. Spurious singularities refer to kinematic regions where the full 
amplitude is nonsingular, but different components of it can contain (cancelling) divergent 
behavior. If we use the unitarity-based method to obtain coefficients of complete loop 
integrals (including associated rational pieces), then the cut terms should be free of spurious 
singularities. If we extract pure (poly)logarithmic expressions, then we must generally add 
rational terms to cancel such singularities. The result of this procedure is referred to as the 
completed-cut term Cn- 

Next, we have to compute a set of recursive rational terms i?^, corresponding to all dia- 
grams in which the shifted legs are attached to different amplitudes. For each arrangement of 
legs, we must sum over the different (complex) factorizations in that channel, schematically 
shown in fig. 1: tree times loop, loop times tree, or tree times tree times a factorization func- 
tion [59] . The factorization-function contribution — which is equivalent to a propagator or 
vacuum-polarization correction in the A/" = case — does not appear for MHV amplitudes, 
and was therefore unnecessary in ref. [27]. 

11 



FIG. 1: Schematic representation of recursive contributions. The labels 'T' and 'L' refer to tree 
and loop vertices. The multi-particle factorization- function contribution (c) does not appear for 
MHV amplitudes. 

Finally, we compute the residues of the rational part of the completed-cut terms, denoted 
by CRn, in the channels affected by the shift. The computation of the residues of CRn{z)/z 
on the physical poles gives us the overlap terms 0„, which correct for double-counting of 
terms between the recursive diagrams and the completed cut. 

The amplitude is the sum of these three terms. 



An = cr 



D 



Cn + Rn + On 



(3.2) 



(See section V of ref. [27] for a relatively simple example of overlap contributions for n = 5.) 
In the present paper we will modify this construction somewhat to allow also for non-trivial 
contributions from z ^ oo. Note that in the present paper, as in the derivation in section 3 
of ref. [27], these individual contributions are defined with respect to A, whereas in the 
explicit calculations in ref. [27], these quantities were defined with respect to pure- finite 
terms {F^ parts of amplitudes). This means the explicitly computed quantities in ref. [27] 
differ from the quantities in the present paper by a factor of i. 

B. Choice of Shifts 

What shift should we choose? The computation in ref. [27] corresponds to choosing a 
[4,5) shift (j = 4, / = 5, in eq. (3.1)) here. However, several properties of the amplitude 
under this shift do not generalize to higher-point amplitudes. In particular, the amplitude 
may not vanish as the shift parameter z is taken to infinity. A [— , +) shift, such as [3,4) or 
[1, 5) in the present case, appears quite generally to have good behavior at infinity. For the 
five-point case, we can verify this explicitly using the known answer [48] , given in eq. (A24) 
of the appendix. For reasons we shall comment on in section IV, for our purposes, here it is 
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FIG. 2: The recursive diagrams arising from a [3,4) shift in A^.^ ^{1 ,2 ,3 ,4^,5"''). Diagram 
(b) has a non-standard complex singularity. 



convenient to introduce a modified L2 function 

Li(r) + I 



UJr) 



1 — r 



U(r) - — 



1r ' 



(3.3) 



ratlier than the more standard L2(r) defined in eq. (A17). These functions differ only in 
terms that are nonsingular as r ^ 1. Using L2a(^), we can give an alternate expression for 
A^.^*^ ^ instead of the form in eq. (A24). 



\M=Q 



/l^r(l-,2-,3-,4+,5 



\Ai=\\-, 2-, 3-, 4+, 5+) + ^crAf ^(1-, 2", 3", 4+, 5+) + cAa 



+ i 



^^ (12) [24] [5 2] (23) ( [5 1] (12) [24] + [5 2] (23) [34] j L^J ffff 



[1 2] [2 3] 



(3.4) 



'51 



where 



R^a = i 



1 (13)'(si2 + S23) ^1(12)(2 3)[2 4][2 5](3-|(2-1)|5-) 



_3 [12] [2 3] (3 4) (4 5) (51) 6 
Consider now the [3,4) shift, 

A3 — > Aq — ZAa , 



[12][2 3](3 4)s2^ 



A4 — ^ A4 + ZA2 



(3.5) 



(3.6) 



The three non-vanishing recursive diagrams are shown in fig. 2. These diagrams correspond 
to residues of the shifted amplitude at poles in z where intermediate states go on shell. 
Diagrams 2(a) and 2(c) are straightforward to evaluate, because the three-point vertex is 
one which appears at tree level, and which can have only a single pole. Diagram 2(b), 
however, involves a one-loop "vertex" A^ (2~ ,3~; K'^). From refs. [30, 31], we know that 
the related "vertex", with opposite intermediate helicity, A^ (2^, 3^; K^), does not factorize 
in complex momenta as a naive generalization of the factorization in real momenta. This 
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FIG. 3: The recursive diagrams arising from a [1, 2) shift in A^.^ ^{1 ,2 ,3 ,4+, 5^). As discussed 
in the text, only diagram (d) is nonvanishing. 

property is related to the appearance of double poles at the loop level. In that case it 
was possible to deduce the relatively simple nonfactorizing structure, at least for the finite 
one-loop helicity amplitudes studied in refs. [30, 31]. For the case of A^ (2^ ,3^; K'^), 
however, we do not know the general structure. Analysis of the behavior under shifting 
of A5^f°(l~, 2~, 3~,4"'", 5^) (see eq. (3.21) below), and of other known amplitudes, reveals 
that it is more subtle than the case of A^ (2~, 3~; K''). (It may even be that, in situations 
where double poles can appear, additional contributions arise which cannot be interpreted 
as factorized diagrams at all. However, an analysis of the diagrams such as those in fig. 2, 
which incorporates some empirical information about the non-standard terms, appears to 
cover any such additional contributions as well.) 

Can we avoid diagrams like 2(b)? To study this, let us consider a [1,2) shift, 

Ai ^ Ai - zX2 , A2 ^ A2 + zXi , (3.7) 

which does in fact avoid generating diagrams whose complex factorization is as- yet unknown. 
The recursive diagrams for this shift are shown in fig. 3. (In the five-point case, choosing a 
[4, 5) or [5, 4) shift would avoid non-standard complex singularities, but as noted above its 
properties do not generalize simply to higher-point amplitudes.) 

Before proceeding to inspect the specific diagrams in fig. 3, we make a few general remarks 
about the properties of three-point vertices, at one loop and beyond, which will be relevant 
for diagrams fig. 3(a) and fig. 3(e). Prior to assigning definite helicities to a three-point 
vertex with external legs a and b, it can be written as A'j^{ea,Sb), where ea,£b are the 
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external polarization vectors, and /i is the Lorentz index for the intermediate gluon. This 
gluon is going on shell in a particular way; either (a h) or [a h] is vanishing, depending on the 
choice of shift. Due to Bose symmetry, and the antisymmetry of the extracted color factor, 
Ag (£„,£;,) is antisymmetric under the exchange ka ^^ h, Sa '^ ^h- Using Bose symmetry 
and gauge invariance, there are only two possible terms in the tensor decomposition of 
Al{ea,e^) [60,61], 

^3 \^ai ^bj = QiySabi (ka+k^)-i)' ^^« ^b ' ^a — £5 5a ■ K^ + /Cj, Ea ' £b) 

Safe 

+ ..K.p;^)*;^(.....-^-^^). (3.8) 

where the form factors gi and g2 are symmetric under ka ^^ ki,. (The required antisymmetry 
in ka ^^ kb follows from the subleading nature of terms proportional to K^ = —{ka + /cfe)'^.) 
We have introduced a fixed external vector rj to indicate that gi and g2 may depend on how 
the intermediate gluon is going on shell. For example, in a real collinear limit, .^ °^'^. is 
the longitudinal momentum fraction carried by gluon a. The form factors can also depend 
on the vanishing quantity Sab- However, the leading dependence can only be logarithmic, 
and so it is subdominant to the power-law behavior of the tensor structures. 
The first tensor structure in eq. (3.8) is the one that appears at tree level, 

^r'^'^l^a'^b) = i^a^h ■ ka - e',^ Ea ' kb + k'^ Ea " Sb) , (3.9) 

^ab 

SO we know a lot about its behavior in complex on-shell kinematics. The second tensor 
structure vanishes for opposite-helicity gluons; with reference vectors qa and g;,, 
^+ ^- e^-kbE^ -kg _ [aqb]{hqa) [ah]{hqa) {ha)[aqb] _ ^ 
" ^ ka-h {aqcDlhqb] {a qa) [h qi] {ah) [h a] 

Therefore, in the case that the two external gluons have opposite helicity, if the tree-level 
vertex vanishes, the loop-level vertex (at any number of loops) should also vanish, since the 
same tensor structure is all that enters. 

If the gluons have the same helicity, say both positive, the second tensor structure is 
nonvanishing off shell, 

5+ ■kbE'l ■ ka [a h] {qb qa) [a h] {h qa) [b a] {a qb) 



el ■ E. 



ka ■ kb {a qa) {b qb) {a qa) {b qb) {a b) [b a 

[ab] 



{ab) {aqa) {bqb) 

[ab] 

{at^' 
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{{ab) {qaqb) + {bqa) {aqt)) 

(3.11) 



However, it vanishes if we approach the complex on-shell kinematics such that [ah] — > 
0. Similarly, the structure relevant when both gluons have negative helicity vanishes as 
(a h) -^ 0. These configurations are those for which the corresponding tree-level vertices are 
also known to vanish. In summary, whenever a tree-level three-point vertex vanishes, the 
corresponding loop-level vertex should vanish as well. 

Note that the identical-helicity tensor structure (3.11), which has a vanishing form factor 
g2 at tree level, but not at one loop and beyond, has the form of an "unreal pole" [31]; that 
is, it is nonsingular for real collinear limits, but blows up or vanishes in complex on-shell 
kinematics. In the case that it blows up, the additional factor of 1/ Sah can produce a double 
pole in (a 6), for the appropriate helicity of the intermediate gluon P. (The contraction 
Sp ■ ka is proportional to either [ab] or (ab), depending on the helicity of P; the former case 
leads to the double pole.) The subleading terms in the expansion around such a double pole 
are non-standard, and are not yet understood. But again, this problem can occur only for 
the identical-helicity case, and only for the complex kinematical configuration for which the 
tree vertex is nonvanishing. 

Now we return to the diagrams of fig. 3. First note that diagram 3(c) contains a tree- level 
vertex which vanishes in the complex on-shell kinematics, because 

Ar%2-,3-,-Kt,)^ (2 3)^ = 0. (3.12) 

As a result, diagram 3(c) vanishes. Note that A^r^^'^lQ,' ,3~ , —K23) would not have vanished 
with a shift of A2 instead of A2. From the above discussion, the loop vertex for the same he- 
licity configuration and type of shift should also vanish. This vertex appears in diagram 3(a). 
Thus diagram 3(a) also vanishes. In the particular case at hand, we can easily verify that 
diagram 3(a) must vanish: If we examine the amplitude (3.4), we see that there is no pole 
in the (23) channel. The reason is that only [2 3] products, not (2 3) products, are present 
in the denominators. The [2 3] products are left untouched by the [1,2) shift. Hence they 
cannot give rise to a pole in z corresponding to diagram 3(a), which would be located at 
(2 3) =0,or z = -(2 3) /(I 3). 

Diagram 3(b) vanishes for the same reason as diagram 3(c). Diagram 3(e) contains a 
loop vertex with a configuration for which the corresponding tree vertex (which appears in 
diagram 3(d)) is nonvanishing. However, because it has opposite-helicity external gluons, we 
know that only the "standard" tree-type tensor structure in eq. (3.8) contributes. The one- 
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FIG. 4: The overlap diagrams arising from an auxiliary [3,4) shift in A^-^ ^{1 ,2 ,3 ,4+,5^ 



loop form factor gi for the A/" = case of a scalar in the loop can then be extracted from the 
one-loop splitting amplitude, and it vanishes. (In the notation of ref. [5], r^ (Ap, a^, b^) = 0, 
for either sign of Xp.) In summary, all the recursive diagrams in fig. 3 are under control, 
and only diagram 3(d) is nonvanishing. 

The [1, 2) shift, however, fails to satisfy the second requirement that An{z) ^ as z ^ oo. 
For large z, we find that 



yl^r°(l",2-,3-,4+,5+;z) 



cr(13)'(z(15) + (2 5)) 



0(i 

.z 



(3.13) 



3 (15)' (3 4) (4 5) [12] 

It is easy to see that the following rational function has the same large-z behavior as the 
full amplitude, 

.cr (13)' (2 5) 



Inf A^f°(l-,2-,3-,4+,5^ 

[1,2> S'-^ ^ ' 



(3.14) 



3 (15)' (3 4) (4 5) [12]' 

where Inf [1^2) ^ denotes those terms in A that would give rise to a non-vanishing contribution 
upon performing a [1, 2) shift and taking z -^ 00. (See section V A for further discussion of 
such terms.) 

We note in contrast that A-^^^ does vanish in the limit. More generally, supersymmetric 
loop amplitudes appear to vanish at large z whenever the corresponding tree-level amplitudes 
do. 

Based on our discussion, it is not surprising that with more positive-helicity gluons, 
higher-point A/" = amplitudes in general suffer from one of these two problems (or perhaps 
both) under any shift: either a 'bad' channel with non-standard complex singularities, or 
'bad' large-z behavior. 



C. Pairs of Shifts to the Rescue 

To get around these problems, we use a combination of shifts, one particular shift to 
obtain one set of contributions, and a different shift to obtain the remaining contributions. 



17 



In our five-point example, we take the [1, 2) shift as the primary one for determining all 
terms except for the large-z contributions given in eq. (3.14). We then use the [3,4) shift 
as an auxiliary shift for determining the large-z contributions of the [1, 2) shift. To see how 
this might be possible, imagine constructing the amplitude ^^^"'^(1", 2~, 3~,4"'", 5"'') via the 
auxiliary [3, 4) shift 

A3 ^ A3 — u'A4 , A4 ^ A4 + w A3 , (3.15) 

where we have taken the shift parameter to be w, to distinguish it from the primary shift 
parameter in eq. (3.7). As above, we decompose the result into three pieces: the completed- 
cut terms, the overlap diagrams, and the recursive diagrams. Let us examine how each 
of these [3,4)-decomposed pieces behaves under the action of the primary [1,2) shift, in 
particular as the parameter z of that shift goes to infinity. To the recursive diagrams for 
the [3,4) shift shown in fig. 2, we add the overlap diagrams in fig. 4. 

In computing the overlap diagrams in fig. 4, we choose a completed-cut expression based 
on eq. (3.4), 

a5(l-,2-,3^,4+,5+) = 

^A^=\l-, 2-, 3-, 4+, 5+) + ^^r (1-, 2-, 3-, 4+, 5+) 
3cr 9 

,: (12) |2 41 |5 2| (2 3) (|5 1| (1 2) |2 41 + [52] (2 3) |3 4|) hj^) 

where we have kept the 2/9^5^''^'^ term in the completed-cut term for consistency with our 
later computations. With this choice, it is a simple matter to check that under the [1,2) 
shift C5(z) vanishes for z ^ oo, using the J\f = 1 amplitude given in eq. (A23), as well as 
eq. (3.3) for L2a(r). 

Now consider the overlap contributions. These are obtained by extracting the residues of 
CR^{w)/w, on the physical poles, 

^(a) = _M, «;(^) = 3, (3.17) 

(35)' [24]' ^ ' 

after applying the [3,4) shift. (The rational part of the completed-cut terms, C-R5, is 
obtained by setting all logarithms to zero in eq. (3.16).) The absence of a (4 5) product in 
the denominator of eq. (3.16) tells us that the first overlap contribution, fig. 4(a), vanishes. 
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The other overlap contribution is straightforward to evaluate, 

O(^) = -(^ + ^)Ar(l-,2-,3-,4+,5+) (3.18) 



, (1 2) [2 4] [5 2] (2 3) ( [5 1] (1 2) [2 4] + [5 2] (2 3) [3 4] j 3,^^ _ ,^^ 



3 [12] [2 3] 2(s5i - S34)2sii ■ 

Applying a [1,2) shift to this expression and taking z ^ 00, we see that this expression also 
vanishes. This is, of course, not surprising, given that it is extracted from the completed-cut 
piece C5, which vanishes at large z. Thus, there are no [1,2) large- 2; contributions arising 
from cut or overlap terms. 

This leaves us with the recursive diagrams to inspect. The first recursive diagram in 
fig. 2(a) is simple to evaluate, yielding, 

B<..=i(l + 5)_^f^^ (3,19) 

V3e 9/ [12] [15] [2 3] [3 4] ^ ' 

The [1, 2) shift of D'^^^ clearly vanishes as 2; ^ 00. 

While we do not know how to evaluate the second recursive diagram, we can nonetheless 
extract its value by 'reverse engineering' from the final answer. That is, we start from the 
final answer in eq. (3.4) (or equivalently, the form in eq. (A24)), shift it, and extract the 
residue of —A-^^^{w)/w at w = [2 3] / [2 4], with the logarithms set to zero. We find the 
contributions corresponding to the sum of diagrams 2(a) and 2(b). Subtracting off the value 
of diagram 2(a), eq. (3.19), gives us the reverse-engineered value for the non-standard residue 
encoded in diagram 2(b), 

^(b) ^ ^ [ (12) (2 3) [2 4] [2 5] (^15 + g34)((5+| 12 |4-) + (5+| 23 |4-)) 
6 1 (15)(3 4)[12][15][2 3][3 4](si5-S34)2 

(12) (2 3) [2 4] [2 5] [4 5] ^ (1 3)^ (1 2) [24] ] 

[12][2 3]s34Si5 [1 2] [2 3] (3 4) (51)' J ■ ^ ' ' 

Using some spinor-product identities, this result can be simplified to 

The appearance of the unreal pole (2 3) / [2 3] in this diagrammatic contribution is the hall- 
mark of a non-standard factorization. If we now shift this expression under the primary [1, 2) 
shift, and take the large-parameter limit, we see that both diagrams 2(a) and (b) vanish. In 
summary, the only piece of A^^^(l~ , 2~, 3", 4+, 5+), decomposed according to the auxiliary 
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[3,4) shift, which survives in the large- 2; hmit of the [1,2) shift, is that given by the [3,4) 
recursive diagram in fig. 2(c). 

Therefore, while the [3,4) shift is not useful for evaluating the entire amplitude, because 
of the non-standard factorization in diagram 2(b), it is useful for deriving a recursion for 
those terms with bad large-z behavior under a different shift. And once we have those terms, 
we can make use of the primary [1, 2) shift to compute the entire amplitude. 

Retaining only those terms from diagram 2(c) that contribute in the large- z limit of the 
[1, 2) shift, we find a remarkably simple recursion formula for these terms, 

Inf <r0(l-,2-,3-,4+,5+) = IniA^=\r,2-,3-,kt,) x — x Ar%-K^„l^,5^) ■ 

11:2) [1,2) S45 

(3.22) 
Here, we denote the operation of extracting the laige-z behavior of the [1, 2) shift by Inf [1^2) • 
We will give a formal definition of this operation in section VA below. In this recursion 
relation, a denotes a shifted momentum with the shift parameter frozen to the value, 

^(c) ^ _M , (3.23) 

(3 5)' ^ ' 

according to the auxiliary [3,4) shift. The right-hand side of the recursion relation (3.22) 
simplifies, because both legs 1 and 2 that are shifted under the primary shift lie on one side 
of the pole, so only the first factor is affected by the extraction of the large- 2; behavior. 

To evaluate the recursion relation (3.22), we need the source of the large-z behavior in 
the four-point amplitude, A^^'^{1^,2^,3^, K^^), which may be obtained by relabeling the 
parity conjugate of eq. (A14) in the appendix. We find, 

Inf<-(l-,2-,3-,4-,5-)^-.^^ ^''^^''^' ' f^^]' 



[1'2> °' 3 (ir45l)[12][2 3] (3i^45) ^45 [i^45 4] [5i^45] 

= ,^ (^3)^(2 5) 

3 (15)' (3 4) (4 5) [12] ' ^ ' ^ 

an expression identical to eq. (3.14). We see that the recursion relation (3.22), containing 
only the (45) channel, reproduces the entire large-z behavior of the five-point amplitude. 

With this term in hand, we can compute the five-point amplitude using the [1, 2) 
shift, even though the amplitude does not vanish as z -^ 00, because the difference 
^'5-r'^ — Inf [1,2) A'5^r^ does have good behavior at infinity. The difference is given by the 
sum (3.2) of the completed-cut expression (3.16), the recursive diagrams (only diagram 3(d) 
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is nonvanishing) , and the overlap diagrams (only one is nonvanishing) . This gives us a 
modified formula for the amplitude, 



An = Inf An + Cr 



Cn + R„ + O 



(3.25) 



where the new term, Inf An, accounts for a non- vanishing contribution at z ^ oo, compared 
to eq. (3.2). (We will need to modify this formula in subsequent sections to account for the 
possibility that the completed-cut terms also have nonvanishing large-z behavior.) 

D. Summary of Strategy 

In summary, the above example suggests a simple strategy for constructing loop ampli- 
tudes, avoiding difficulties from either non- vanishing large- 2; behavior or from non-standard 
complex factorizations: 

• Use a primary shift whose recursion relation does not contain any non-standard com- 
plex factorizations, but under which the amplitude might not vanish at large z. 

• Use an auxiliary shift and recursion relation to determine the laige-z behavior of the 
amplitude under the primary shift. This auxiliary recursion relation may contain non- 
standard complex factorizations, but these will be harmless if the contributions from 
these channels vanish in the large-z limit of the primary shift. 

In general, of course, we may not know ahead of time whether an amplitude has non- 
trivial large-parameter behavior for a chosen shift. We should therefore first derive an 
auxiliary recursion for this behavior (which might of course reveal that it is absent). The 
derivation of an auxiliary recursion relation may require assumptions about the behavior of 
contributions in 'bad' channels, such as that of the vanishing of diagram 2(b) in the large- 
parameter limit of the [1, 2) shift. These assumptions can be verified once we have obtained 
a (candidate) final answer for the full amplitude, through a check of all (real-momentum) 
factorization limits. If the assumptions are not valid, and the auxiliary recursion relation 
yields an incomplete expression for the large-parameter behavior, the final answer will not 
have the correct collinear or multi-particle factorization limits. In section VI, we present 
shifts that we expect will satisfy all criteria necessary for constructing the rational parts 
of any n-gluon amplitude. We also expect this strategy to be widely applicable to other 
phenomenologically interesting amplitudes. 
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IV. A NON-MHV SIX-POINT EXAMPLE: ^6(1-, 2", 3-,4+, 5+,6+) 

In the previous section we have developed a strategy for computing amphtudes even 
with shifts under which An{z) y^ as the shift parameter z ^ oo. We now apply these 
ideas to the calculation of a previously unknown amplitude with three negative helicities, 
<r0(l",2^,3~,4+,5+,6+). 

As in the five-point case, we shall use a [1,2) shift as the primary shift for computing the 
amplitude. For the same reasons as in the five-point case, the six-point amplitude is also 
free of contributions in 'bad' channels with non-standard complex singularities. So we can 
evaluate all of the recursive diagrams. It is clear from a consideration of the factorization 
properties as the momenta of gluons 5 and 6 become collinear that the amplitude must have 
non-trivial behavior at large z, because the corresponding five-point amplitude has such 
behavior under the shift. Our first task is to obtain a rational function which reproduces 
that behavior. 

The full amplitude is obtained by combining the large-z terms with the other standard 
terms — completed-cut terms, recursive diagrams, and overlap diagrams — via eq. (3.25). 
It turns out that we need to modify eq. (3.25) a bit more to subtract out nonvanishing 
large- 2; contributions from the completed-cut terms, C„. As we shall discuss more fully in 
section V, we are really performing a contour-integral analysis on the rational function of 
z, {An{z) — [Inf A„](2;)) — Cn{z). It is this expression that must vanish as 2; ^ 00 to allow 
closing the contour at infinity. In the five-point example of the previous section, we avoided 
the need for subtracting the large- 2; behavior of C„, by use of the modified L2 function in 
eq. (3.3); we simply ensured that C^lz) ^ as z ^ 00. For more general amplitudes it 
is simpler to subtract out any nonvanishing large-z contributions arising from C„, before 
adding back the complete large- 2; behavior. This amounts to performing a contour-integral 
analysis on {An{z) — [Ini An]{z)) — {Cn{z) — [Ini Cn]{z)) , which automatically vanishes at 
infinity. The full result is then. 



An = Inf An + cr 



Cn-lnfCn + Rt^ + 



n 



(4.1) 



where Inf A„ contains the large-z behavior of the entire amplitude, R^ are the contributions 
from the recursive diagrams, and 0„ are the overlap pieces. In section V, we will provide a 
more systematic derivation of this formula. 
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FIG. 5: The recursive diagrams arising from an auxiliary [3, 4) shift in 

Aq[^^{1^ ,2~ ,3~ ,4"^ ,5^ ,6^). Diagram (c) has non-standard complex singularities. 

A. A Recursion Relation for Large-Parameter Behavior 

We begin by conjecturing a recursion relation for the large-z behavior under the [1, 2) shift 
that is the simplest possible generalization of the one we found for the five-point amplitude, 
namely, 

Inf<r0(l-,2-,3-,4+,5+,6+) = Inf <r°(l-> 2", 3", i^+, 6+) x — 

[1>2) [1,2) S45 

x/ir(-K4-5,4+,5+). (4.2) 

As in eq. (3.22) for the five-point case, a indicates a momentum shifted according to the 
auxiliary [3,4) shift in eq. (3.15), with w frozen to the value 

..™ = -M. (43) 

(35) '• 

We will return to a more systematic construction of the large-z contributions via on-shell 
recursion in sections V-VII. Here we only wish to provide some motivation for the large-^ 
recursion relation (4.2), following the same procedure as at five points. 

The recursive diagrams for computing the amplitude using a [3,4) shift are displayed in 
fig. 5. However, we are only interested in the values of the diagrams in the large-z limit of 
the [1, 2) shift. Diagram 5(a) has legs 1 and 2 attached to a tree which vanishes in the large-^ 
limit [29]. As we shall discuss in section VD, we assume, based on empirical evidence from 
known amplitudes, that the full diagram is also suppressed at large z, even though the loop 
vertex is non-standard. Diagram 5(b) is more complicated, in that the shifted legs straddle 
the pole, but its value may be determined easily since its components are known; it is then 
not difficult to verify that it is suppressed at large z. Diagram 5(c) is problematic, since it 
also contains a non-standard complex pole. For this case, we assume that the addition of 
the extra leg on the tree side, compared to the corresponding five-point diagram shown in 
fig. 2(b), does not upset the suppression at large z; we take it to have a vanishing contribution 
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in this limit. One must also check that there are no additional large-z contributions from the 
cut terms (after subtracting Inf C„) or from the overlap terms. Following a similar analysis 
as for the five-point case, it is not difficult to verify that there are none in this case. We are 
left with the single diagram 5(d), as the sole surviving contribution in the large- z limit of 
the primary [1,2) shift, motivating eq. (4.2). 

To evaluate the large-z recursion relation in eq. (4.2), we use the five-point large- 2; result 
in eq. (3.24) as input. We obtain, 

,Ar=on- o- .-.+ .+ .+^ _ .^r (13)=^ (26) 1 [45f 



Inf A^f°(l-,2-,3-,4+,5+,6+) = -2— ,^-^. . ^^^ ^ 

[1.2> '•' ' 3 (16)2 (3i^45) (i^45 6)[12] S45 [i^45 4] [5i^45] 

= ,££ (^3)^(2 6) 

3 (16)' (3 4) (4 5) (5 6) [12]' ^ ' ^ 

where, as above, a indicates a momentum shifted according to the [3,4) shift in eq. (3.15), 
with w frozen to the value (4.3). 

B. The Completed-Cut Terms 

The cut-containing terms of our target amplitude were computed in ref. [26] using a recur- 
sion relation on coefficients of integral functions, along with known lower-point results [6, 48]. 
In the six-point case, this procedure yields, 

56(l-,2-,3-,4+,5+,6+) = ^A6^r'(l">2-,3-,4+,5+,6+) 

+ ^Ar (1-, 2~, 3-, 4+, 5+, 6+) + C^, + C^, , (4.5) 

y flip 1 

where 

1 2) (2 3) [2 4] (1- 1 (3+4) 12") [(3"! 42 |1+) S234 - (3"! 2(3 + 4) |1+) S34] ^2{ 



c: ' 



-'6 



-^234 ' 
-S34 ' 



(34)(56)(61)[23](5-|(3+4)|2-) s^ 

(3 5) [4 5] [5 6] (5- 1 (1 + 2) |6-) [(3-| (5 - 4) |6-)s345 + (S"! (4 + 5) 16^)534] U 



34 
S345 "n 
-S34 



.3 



(4 5) [1 2] [16] (5- 1 (3 + 4) |2-) si, 

(4.6) 



and where we have introduced the fiip symmetry operation. 



X(l,2,3,4,5,6) =X(3,2,1,6,5,4). (4.7) 

flip 1 
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The first term in eq. (4.5) is proportional to the contribution of an A/" = 1 chiral multiplet 
in the loop. This contribution is fully constructible from the four-dimensional cuts [6]. The 
result is [24], 

(4.8) 



A^r\^~, 2-, 3-, 4+, 5+, 6+) = S^, + S^, 



flipl 



where 

icr 






6 



2 



^.4^(1", 2-, 3-, 4+, 5+, 6+) Ko(s34) (4.9) 



t 



(1-|(2 + 3)|4')^ [(3-|42|l+)s234-(3-|2(3 + 4)|l+)s34]Lo(^) 

(5 6) (6 1) [2 3] S234 su (5- 1 (3 + 4) |2-) su 

(3-| (1 + 2) |6-)^ [(3-| (5 - 4) |6-) ^345 + (31 (4 + 5) |6-) su] U{^ 
(3 4) (4 5) [12] [16] S345 (5- 1 (3 + 4) |2-) S34 

After performing the [1, 2) shift, the completed-cut expression given in eq. (4.5) does not 
vanish as 2; — > 00, but tends to a purely-rational constant, 

Inf deir, 2-, 3-, 4+, 5+, 6+) = lim Cair , 2", 3^,4+, 5+, 6+; z) 

[1,2) 2-+00 



z 



(1 2) (1 3) (3^1 (4 + 5) |2-) - (1-| 2(4 + 5) |3+) + (1 3) S345 



6 [1 2] (3 4) (4 5) (6 1)' S345 (5~ | (3 + 4) |2-) 

^ (12) [2 4] (13) ((12) [2 4] -(13) [3 4]) 

^6 [2 3](5 6)(61)S34 (5-|(3 + 4)|2-) ' ^ ' ' 

Since we have already determined the complete large-z behavior of the full amplitude in 
eq. (4.4), we must subtract from Cq this rational constant, Inf [1^2) Cq, so that the difference 
Cq — Inf [1^2) C'e vanishes as z —^ 00. When computing the overlap contribution we may use 
either this difference or the original Cq. They are equivalent because Inf[i^2) C*6 has no poles 
in z under a [1,2) shift; therefore it does not generate an overlap contribution. 

C. Recursive Contributions 

Next we evaluate the recursive diagrams for the [1, 2) shift oi A^^'^(l^ , 2~, 3~, 4+, 5+, 6+). 
Most of them vanish. Some of the vanishing diagrams are shown in fig. 6. (We omitted two 
diagrams in the (23) channel, where (— -^^23) carries negative helicity, which vanish even 
more trivially.) The first two diagrams, 6(a) and 6(b), vanish because the loop vertices 
i?3(6"^, 1^, — A'^) vanish, as discussed in section IIIB and in the appendix. Diagrams 6(c) 
and 6(d) vanish because 

Af^(6+,i-,-A:+)oc [6A:6i]' = 0, A'r{2-,3-,-K+,)^ (2 3)^ = 0. (4.11) 
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FIG. 6: Some vanishing recursive diagrams for the [1, 2) shift of ^gf]^ '^(1 ,2 ,3 ,4+, 5^,6+). 




FIG. 7: Non-vanishing recursive diagrams. Diagram (i) is the factorization-function contribution. 

As discussed for the five-point case, diagram 6(e) vanishes because its loop vertex is of the 
same type as the vanishing tree vertex in diagram 6(d). Summarizing, we have 



D^ 



(a) 



D, 



(b) 






Dlt'^ =0. 



The four non- vanishing recursive diagrams, 



(4.12) 



(4.13) 



are shown in fig. 7. These diagrams are straightforward to evaluate because all channels 
involve standard factorizations. Diagram 7(f) yields 



Af^il-, ~K^„ 6+) — R,{2', 3-, 4+, 5+, K, 



61V 



-.r^+^v 



(3-|(l + 2)|6-)^ 



3e 9; [12] (3 4) (4 5) [61] S345 (5-| (3 + 4) |2-) 

^ ^345 + ^34 (3 5) [4 5] [5 6] (5- 1(1 + 2) |6-) 

6 534^345(5345 - 534)2 [1 2] (4 5) [6 1] (5- 1 (3 + 4) 12") 
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x((3-| 4 |5-) (5-| (1 + 2) |6-) + (3-| 5 |6-) S345) 
i [4 6]^(4-|(3 + 5)|2^) ^ i [46]^ (3-| (1 + 2) |6-) 



3 [12] [2 3] [3 4] (4 5) [61] (S^l (3 + 4) |2") 3 [1 2] [3 4] (4 5) [6 1] (S"] (3 + 4) 12") 
i (35) [45] [56] (3^1 (1 + 2) 16") (5^1 (1 + 2) 16") 



6 [12](4 5)[61]S34S345 (5-|(3 + 4)|2-) 

where the five-point scalar loop vertex is given in eq. (A26). Diagram 7(g) is 



(4.14) 



z 



Dfe) = Ar(2-,3-,4+,i^+J i?4(i-,-i^56i>5+,6+) 

S561 

= (^ + H)^r(2-, 3-, 4+, i^+,) — ^r (i-, -Kei, 5^ 6+) 
V3e 9/ S561 

'l + ^^ (1-1(2 + 3)14-)^ (4,5) 

.3e 9/ [2 3] [3 4] (5 6) (61) S234 (5-| (3 + 4) 12") ' ^ ^ 

where we used the four-point tree amplitude, eq. (A4), and the loop four- vertex given in 

eq. (A21). It is easy to see that diagram 7(h) gives the same value, 

^(h) _ ^(g) ^ (4_,g) 

Diagram 7(i) contains the factorization function contribution [59], which for the scalar 
loop case amounts to a vacuum polarization insertion. The value of the factorization function 
vertex appearing in this diagram is given in eq. (A12). Using this value of the factorization 
function, diagram 7(i) is given by 

V3e 9/ S561 

= -Dfe). (4.17) 

The fact that diagrams 7(g), (h) and (i) are equal, up to signs, is rather special to this 
amplitude. 

D. The Overlap Contributions 

To evaluate the overlap contributions we start from the rational parts of the completed- 
cut terms, obtained by setting all logarithms to zero in eq. (4.5). That is, we replace the Kg 
and Lj functions defined in eq. (A17) with their rational parts, 

1 



Ko(r) 



e 



Lo(r) ^ 0, 

1 (r + 1) 
Uir) -^ - , • (4.18) 

A J 2r 1-r 2 ^ ' 
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FIG. 8: The overlap diagrams for the [1,2) shift. Diagram (a) vanishes. 



Making these replacements in eq. (4.5) for Cg gives us CRq. Applying the [1,2) shift, 
eq. (3.7), yields CRq{z), which we use to evaluate the overlap contributions. 

As depicted in fig. 8, for the [1,2) shift there are three channels which can potentially 
contribute to the overlap, 

O = O^'^) + O^^) + O^'^) , (4.19) 

corresponding to the three residues of CRq{z)/z located at the values of z, 

(2 3) .b) [16] rc^ S234 



.(a) 



.(c) 



(13)' ^ [2 6]' ^ (1-|(3 + 4)|2-) • ^^'^^^ 

Evaluating these residues gives us the overlap contributions. After simplification, they are 
given by, 

0(a) 







0(b) 



h. 



(3- 



'1 



2) |6-)' 



3e 97 [12] (3 4) (4 5) [61] 5345 (5-| (3 + 4) |2-) 
i (1 2) (3^1 (4 + 5) |2~) (3-| (1 + 2) |6-)' 
'6 [12](3 4)(4 5)s6iS345 (5-| (3 + 4) 12") 

t_ [4 5] [5 6] (3 5) (5- 1 (1 + 2) |6-) ((3 5) [5 6] S345 - (3 4) [4 5] (5"! (3 + 4) 16")) 
6 



[1 2] (4 5) [61] (5- 1 (3 + 4) |2-) 



X- 



-534 + "5345 



5345345 (5345 '~ S34) 

i (15)[4 5][4 6](1-|(5-6)|4-) 
'6[2 3][3 4](5 6)s6i(5-|(3 + 4)|2-) 



0(c) 



—t 



I 



h. 



(l-|(2 + 3)|4-)= 



3e 9 J [2 3] [3 4] (5 6) (6 1) S234 (5" | (3 + 4) |2-) 



-\2 



+ 



+ 



(12)[24](1-|(2 + 3)|4-) 



6 [2 3][3 4](5 6)(61)S234 (5-| (3 + 4) 12") 
2 (15)(34)[45](16)(1-|(2 + 3)|4-)' 



(4.21) 



6 [2 3](5 6)(61)=^S34S234 (5-|(3 + 4)|2-)' 
Although there is no need to keep the 1/e terms (their values are known a priori), we have 
carried them along. 
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E. The Full Amplitude and Consistency Checks 



We may now combine all the pieces to obtain the full amplitude, 



\M=0/ 



Cr + Rr 



(4.22) 



A^e7^(l-,2-,3-,4+,5^6+) = cr 

where Cq is given in eq. (4.5). The rational remainder i?6, consisting of recursive diagrams, 
overlap diagrams and large- 2; contributions, is, 

^6 = — Inf A^r° - Inf Ce + /^^^^ + ^^^^ + ^^^^ + ^^'^ + O(^) + 0(=) , (4.23) 

Cr [1,2) ^'' [1,2> 

where Inf[i^2> ^ea ° i^ given in eq. (4.4) and Inf^i 2) Cq in eq. (4.10). The values of the recursive 
and overlap contributions are given in eqs. (4.14)-(4.17), and eq. (4.21). Simplifying the 
result for the rational remainder and making use of the flip symmetry (4.7), we can write 
the result for the remaining rational part in a rather compact form. 



D DO I oa 



flipl 



where 

i 1 r [4 6]^ [2 5] (5 6) (13)^ (2 5) [2 3] 



(4.24) 



^^ 6[2 3](5 6) (5-|(3 + 4)|2-) i [12] [3 4] [61] (3 4) (4 5) (61) 

(l-|(2 + 3)|4-)V(l'|2|4-)-(l-|5|4-) (13) [46]x 

[34] (61) I S234 (34) [61]/ ^ • ^ 

(1 3)^(3 (1~| 2 |4-) + (1^1 3 |4-)) [46]^(3 (l^l 5 14") + (1~| 6 |4-)) 

The result (4.24) is manifestly symmetric, not only under the flip (4.7), but also under the 
second flip symmetry, involving spinor conjugation, 

X(l, 2, 3, 4, 5, 6) = X*(6, 5, 4, 3, 2, 1). (4.26) 

flip 2 

The remarkable simplicity of the rational remainder is rather striking. 

We have performed a number of checks on our result for the amplitude, eq. (4.22). We 
have conflrmed that it has the proper factorization properties in real momenta in all two- 
and three-particle channels and that all spurious singularities indeed cancel as they should. 
Since the large- 2; contribution in eq. (4.4) contains kinematic poles in a variety of channels, 
an omitted piece would necessarily be detected in some of the collinear limits. Finally, the 
numerical value at one phase-space point agrees with that in ref. [3]. The consistency of the 
amplitude demonstrates the validity of our procedure for determining the large-z terms, for 
a new and rather non-trivial analytic amplitude. 

29 



V. ON-SHELL RECURSION RELATIONS FOR LOOP AMPLITUDES 

Before continuing to the case of general helicities, we present in this section a more 
systematic discussion of loop-level on-shell recursion relations, emphasizing the extensions 
beyond ref. [27]. The derivation of such loop recursion relations is similar in spirit to the tree- 
level case, but it does require the treatment of factorizations which differ from the 'ordinary' 
factorization in real momenta. It also differs because of the appearance of branch cuts and 
spurious singularities associated with logarithms or polylogarithms. The cut-containing 
parts of the amplitude are an input to the loop recursion relations. We assume that they 
have been computed by other means, such as the unitarity-based method. 

A. Analytic Behavior of Shifted Loop Amplitudes 

The starting point for the loop recursion relations is a complex-valued shift of the mo- 
menta of a pair of external particles in an n-point amplitude, kj -^ kj{z), ki — > ki{z). We 
describe a [j, I) shift in terms of the spinor variables A and A defined in eq. (2.3) , 

[j, I) : \j -^ \j - z\i , Xi ^ \i + zXj , (5.1) 

following the tree-level construction [29]. The shift maintains overall momentum conser- 
vation as well as the masslessness of the external momenta, k'j = kf = 0. Let us denote 
the original n-point amplitude by An = ^ra(O), and the shifted one by An{z). We seek an 
equation for A„(0) relying on the analytic properties of An{z). (We will also denote by f{z) 
other functions / of the momenta, such as the cut-containing terms, after the shift (5.1).) 

At tree level, An{z) is a meromorphic function of z. Its poles are determined by the 
factorization properties of An as multi-particle invariants or spinor products vanish. The 
former are identical for real and complex momenta, and at tree level, the singularities for 
spinor products of complex momenta are also completely determined by the corresponding 
singularities (collinear factorizations) for real momenta. At tree level, one can show [29, 44, 
62], that there are always choices of shift momenta for which An{z) ^ as z ^ oo. This 
allows the derivation of a recursion relation through consideration of a contour integral on 
a circle at infinity. 

At one loop, we must consider several new aspects. The most obvious of these is the 
presence of branch cuts, which arise from logarithms or polylogarithms in the amplitudes. 
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But there are a number of other important features. While factorization in real momenta 
is understood [5, 59, 61], this does not completely determine the singularity structure in 
complex momenta. In particular, we must in general handle double poles as well as 'unreal' 
poles, present for complex momenta but absent for real ones. We obtained some of the 
required factorizations heuristically [30, 31], confirming them by explicit calculation. In 
certain two-particle channels, however, the structure of complex factorization is not yet 
completely clear; accordingly, we will design our calculations to avoid relying on them. 

While it does appear in general possible to find momentum shifts under which An{z) -^ 
as z -^ oo, it is not always possible to do so while avoiding poles in channels with obscure 
complex-factorization behavior. Accordingly, we must deal with amplitudes that have either 
a pole at infinity, so that An{z) ^ oo as 2; ^ oo, or that behave as a finite constant. We will, 
however, assume that this behavior is given by a rational function of z. In all cases we deal 
with here, where the corresponding tree amplitudes vanish in the large-z limit, the large-^ 
terms are purely rational. For the moment, let us assume we know these terms, along with 
the cut-containing terms; below we return to the issue of their computation. 

More precisely, we wish to define an operator Inf that yields a pole-free rational function 
reproducing the laige-z behavior of an amplitude, when both are shifted. 



lim 

2— »oo 



Inf A 



z)-Aniz))=0. (5.2) 



We implement this operator via a series expansion of the amplitude around the point z = 00, 
or u = 1/z = 0. We take An{z) to be an analytic function with a series expansion of the 
form, 

'max ^ 00 00 00 

An{z) = Y,- °^"^ +Y.^' ^^'^ + Yl ^^(")"' ^'^ + Yl ln'(M)M* S'^ , (5.3) 



i=0 i=l i=l i=l 



where the coefficients, a^~'^\ b^^\ c*^*-* and d^^^ depend on spinor products. The leading degree 
of large- 2; behavior, imax, depends on the helicity configuration. In writing this expansion, 
we have assumed that no ln{z) factors appear in terms which survive as z ^ 00. (We 
can always confirm the validity of this assumption for a given amplitude, since all logs and 
polylogarithms will have already been computed via the unitarity method or other means.) 
The large-z terms are those which do not vanish as 2; — > 00, 

A^^^ '{z)=J2 ^' «^"'^ = [ Inf An] (z) , (5.4) 



j=0 
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FIG. 9: A configuration of poles and branch cuts for a term in a one-loop amplitude, with a 
branch-cut-hugging contour. 

given by shifting their progenitor Inf A^. 

We may then write An{z) as a sum of the large-z terms and the remaining terms, 

A^iz) = A'^riz) + Al°'^^iz), (5.5) 

so that all finite- 2; poles, logarithms, and polylogarithms in z are put in the second term. 
We will call the latter object the 'finite-pole' terms. 

Since we will be interested in the unshifted physical amplitude obtained by setting z = 0, 
we will need those terms in A„(0) that generate A^^^ ^{z) under the shift. We find in practice 
that these are given by setting z = in A^^^ ^{z), that is keeping only the 2;° = m^ terms in 
the expansion of An{z), 

lniAn = a^^K (5.6) 

A key aim of this paper is to provide a means for computing this quantity. 

Since y4P°^'^^(2;) vanishes for z — > 00, we may apply Cauchy's residue theorem. Consider 
the contour going along the circle at infinity, but avoiding the branch cuts by integrating 
inwards along one side, and then outwards along the other, as shown in fig. 9. We route the 
branch cuts so that no two overlap. (As discussed in ref. [27] the case of poles touching the 
ends of branch cuts works in the same way as the basic case, and we shall not distinguish it 
further in our discussion.) The integral along this contour is given by the sum of residues. 
This contour, however, includes a branch-cut-hugging integral. 



f ^Ar'-(.) + -^/ -Ar>-(.), (5.7) 
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where B^ is directed from an endpoint Bq to infinity, and B^ is directed in the opposite 
way. Now, A^°^'^^{z) has a branch cut along B, which means that it has a non- vanishing 
discontinuity, 

27ri DiscB AP^^^^z) = Ap°''=^(z + ie) - Ap°1''(^ -ie), z on B . (5.8) 

(Because Aj^§ ^ is taken to be purely rational, Disc^ ^^"^'^^(-z) = Disc^ A„(2;).) Thus, using 
the vanishing of A^°^'^^{z) for z -^ oo, we have, 

E/lpolcs/ \ roo T 

Res " ^^ + / — Discs ^r'"(-2) • (5-9) 

, ^=^« Z Jbo Z 

poles OL " 

B. Cuts and Cut Completion 

To proceed further, let us assume that we have already computed all terms having branch 
cuts, along with certain closely related terms that can generally be obtained from the same 
computation. That is, we have computed all polylog terms, all log terms, and all tt^ terms. 
As discussed in ref. [27], there are also certain classes of rational terms that are natural to 
include with the cut-containing terms. 

In particular, there are rational terms whose presence is required to cancel spurious 
singularities in the (poly) logarithmic terms. Such spurious singularities arise in the course 
of integral reductions. They are not singularities of the final amplitude, because they are 
unphysical, and not singularities of any Feynman diagram. A simple example comes from a 
'two-mass' triangle integral for which two of the three external legs are off shell (massive), 
with momentum invariants s\ and S2, say. When there are sufficiently many loop momenta 
inserted in the numerator of this integral, it gives rise to functions such as, 

where r is a ratio of momentum invariants (here r = Si/s2)- The limit r ^ 1 (that is, 
Si -^ s-i) is a spurious singularity; it does not correspond to any physical factorization. 
Indeed, this function always appears in the amplitude together with appropriate rational 

pieces, 

ln(r) + 1 - r 

(1-r)^ ' ('•'') 

in a combination which is finite as r — > 1. From a practical point of view, it is most convenient 

to 'complete' the unitarity-derived answer for the cuts by replacing functions like eq. (5.10) 
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with non-singular combinations like eq. (5.11). Such completions are of course not unique; 
one could add additional rational terms free of spurious singularities. 

The reason we want to eliminate spurious singularities from sub-expressions has to do 
with the sum over pole residues in eq. (5.9). The sum runs over all poles, whether they arise 
from a shift in a physical singularity variable, or in a spurious one. In practice, it is sufficient 
to eliminate all spurious singularities that acquire a z dependence under the momentum shift. 
By construction, J^^ '''{z) is free of such 'dangerous' spurious singularities, since it cannot 
contain poles in z. 

Singularities that look like collinear ones, but involve non-adjacent momenta, are also spu- 
rious. For example, in the scalar contributions to the five-gluon ( — I 1-+) amplitude [48], 

there are factors of (2 4) and (2 5) appearing in the denominators of certain coefficients. 
These might appear to give rise to non- adjacent collinear singularities in complex momenta; 
but by expanding the polylogarithms and logarithms in that limit, one can show that these 
singularities are in fact absent in the full amplitude. 

Let us define two decompositions of the amplitude. The first is into 'pure-cut' and 
'rational' pieces. The rational parts are defined by setting all logarithms, polylogarithms, 
and TT^ terms to zero, 



Rr,{z\ 



1 

cr 



-/I,. 



-K 



rat 



Cr 



ln,Li,7r2^0 



(5.12) 



(The normalization constant cp, defined in eq. (2.15), plays no essential role in the following 
arguments, we carry it along for completeness.) The 'pure-cut' terms are the remaining 
terms, all of which must contain logarithms, polylogarithms, or tt^ terms. 



CJz\ 



1 
Cr 



- Ar. 



- Ar. 



IrijLijTT^ 



(5.13) 



(5.14) 



pure— cut Cp 

In other words, 

An{z) = Cr Cniz) + Rniz) 

where we have explicitly taken cr outside of Cn{z) and Rn{z). 

The second decomposition uses the 'completed-cut' terms, obtained from Cn{z) by replac- 
ing logarithms and polylogarithms by corresponding functions free of spurious singularities 
(at least those that suffer a shift). We call this completion C„. The decomposition defines 
the remaining rational pieces i?„. 



AJz) 



Cr 



CJz) + RJz) 



(5.15) 
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This reorganization has effectively moved some of the rational terms into the completed-cut 
terms. In general, C„(z) and C„(z) will not vanish as z ^ oo. We will assume that the 
large-z behavior is given by a rational function of the spinor products, as in the previous 
examples. The rational function may include contributions from series expansions of the 
logarithms and polylogarithms in Cn{z). Taking this into account, we can define a useful 
decomposition of the finite-pole terms, 

Al^'^'^z) = cr[Cn{z) + R.^{z) - [\niCn]{z) - [Inf i?„](^)] . (5.16) 

In the cases relevant to the present paper, [Inf C„](-2) is in fact at worst a constant in z. 
Note that the additional terms will not contain any 'dangerous' spurious singularities, as 
the act of taking the large- 2; limit will eliminate them. 

We also need to define the rational part of the completed-cut terms, CRn{z). We write, 

Cn{z)=Cn{z)+CRn{z), (5.17) 

where 

CRr,{z) = dr,{z) . (5.18) 

rat 

Combining eqs. (5.14), (5.15), and (5.17), we see that the full rational part is the sum of the 
rational part of the completed-cut terms, and the remaining rational pieces, 

Rn{z) =^n{z) + Rn{z) . (5.19) 

Now, because we know all the terms containing branch cuts, we could compute the 
branch-cut-hugging integral in eq. (5.9), 

/ - DiscBAP°'^^(^) = / - DiscB[a(^)-[Infa](^)] 

Jbo ^ Jbo ^ 

C^ dz - 

= / — DiscB Cn{z) , (5.20) 

Jbo z 

where the second line follows from the rational nature of [Inf C„](2;). However, there is no 
need to do the integral explicitly, because we already know the answer for the integral, plus 
the associated residues. Up to a contribution coming from Ini Cn, it is just (^^(O), part of 
the final answer. That is, applying the same logic to Cn{z) — [Inf C„](2;) as was applied to 
A^^''''{z) in eq. (5.9), we have, 

C„(0) = Inf C„ - y Res ^^^^^ - / — Disc^ Cn{z) , (5.21) 

poles a ^^0 
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where by construction Inf C„ does not contribute to the sum over residues either. 

Using eq. (5.9), the decomposition (5.15), and eq. (5.21) to evaluate the terms involving 
Cn{z), we can write our desired answer as follows. 



AP°^-(0) = -cr 



oo 



- Disc^a„(^)+ y Res^^+ y Res^^^ 

'- " poles a poles a 



Cr 



a„(0)-lnfa„- y Res^"^""^ 

poles a. 



(5.22) 



By construction, the completed-cut terms Cn{z) contain no spurious singularities, and so the 
sums over the poles in eq. (5.22) are only over the genuine, 'physical' poles in the amplitude. 
As we are working with complex momenta, these are the poles that arise for such momenta, 
and not merely those that arise for real momenta. In particular, this means that double 
poles and unreal poles may appear, as discussed in detail in refs. [30, 31]. 

C. Residues of the Remaining Rational Pieces Rn{z) 

Our next task is to evaluate the residues of the remaining rational terms, -R„. We will 
do this by setting up an on-shell recursion relation. As discussed in ref. [27], pure-cut and 
rational terms factorize independently, and so one can use factorization to construct an 
on-shell recursion relation for the full rational terms /?„. However, there is no fundamental 
factorization distinction between the rational terms in Cn and those in i?„, so we cannot set 
up a direct recursion relation for them. We will instead compute them indirectly, by first 
computing the full rational terms, and then subtracting terms which are present in both the 
full rational terms and in the completed-cut terms. These are exactly terms coming from 
CRni which we therefore call 'overlap' terms, 

- y Res^ = - y Res^^+ V Res ^^ . (5.23) 

^-^ Z=Za Z ^-^ Z = Za Z ^-^ Z = Za Z 

poles a poles a poles a 

Since we know CRn explicitly, we can compute the last sum by shifting and extracting poles, 

O^^y Res ^^^ . (5.24) 

^ ' Z = Za Z 

poles a 

To obtain the first term on the right-hand side of eq. (5.23), we must analyze the poles 
in Rn{z), that is to say its properties at appropriate null complex momenta. The behavior 

36 



of Rn{z) can be extracted from the factorization properties of the amphtude as a whole, by 
following the analysis of ref. [27], and separating two classes of terms — pure-cut and rational 
— in the factorized amplitudes. Only rational terms in the factorization can contribute to 
the required sum of residues. 

Given the [j, /) shift (5.1), we define a partition P to be a set of two or more cyclicly- 
consecutive momentum labels containing j, such that the complementary set P consists of 
two or more cyclicly-consecutive labels containing /: 

P = {P,,P2,...,j,...,P-i}, (5.25) 

P = {Pi, P2,..., /,..., P-i}, 
P U P = {l,2,...,n}. 

This definition ensures that the sum of momenta in each partition is z-dependent, so that 
it can go on shell for a suitable value of z. The complex on-shell momenta kj, ki and P are 
determined by solving the on-shell condition, P^ = 0, for z. 

At one loop, there are in general three contributions to factorization in any given channel, 

A 1-loOp 



^1-loop ^ _J__ ^ ^trce I ^tree ^ ^ ^ ^1-loop ^trec ^ '■ -^ ,, 4 tree 

s ^ S- • R ^ L R ^ L 

In the first two terms, one of the factorized amplitudes is a one-loop amplitude and the other 



is a tree amplitude. The last term contributes only in multi-particle channels, and contains 
a one-loop 'factorization function'. For the case of a scalar in the loop (A/" = 0), this function 
is equal to the scalar contribution to the gluon vacuum polarization [59]. Accordingly, in 
addition to the sum over channels, we will have a sum over these different factorization 
contributions. Taking the rational parts, we obtain, 

Res^^^ = P^(A;i,...,A;„) 



^^ Z=Za Z 

poles a 



J2 $^|i?(A:p„...,^„...,A;p_„-P'^) 

partitions P /i=± ^ 

i 



X X A [kp^ , . . . ,kl, . . . , kp_^ , P ) 



+ A'^'%kp,,...,k,,...,kp_„-P') (5.27) 

^ "p2 ^ -^(^Pi 1 ■ ■ ■ il^h ■ ■ ■ 1 ^"P_i 1 P ) 

+ A''''^{kp„...,k„...,kp_„-P'') 



X 



iRr{P) 



p2 
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xA [kp_^, . . . , ki, . . . , kp__^, P )> 



where P^ is the squared momentum associated with the partition P, and Rjr is defined in 
the appendix. The hatted variables are given as usual by shifting momenta according to 
eq. (5.1), and freezing z to the value that puts P on shell, 

The result (5.27) follows directly from the general factorization behavior of one-loop ampli- 
tudes, plus the separate factorization of pure-cut and rational terms [27]. Although the R 
functions are not complete amplitudes, they can be thought of as vertices. Equation (5.27) 
then gives rise to a set of 'recursive diagrams'. 

Inserting the result into eq. (5.22), and thence into eq. (5.5), gives us the basic on-shell 
recursion relation for complete one-loop amplitudes, 

A„(0) = Inf An + cr p„(0) - Inf Cn + R^ + 0„] . (5.29) 

To compute with this equation, we construct R^ via recursive diagrams; that is, via 
eq. (5.27). The 'overlap' terms 0„ can also be given a diagrammatic interpretation, as- 
sociating each pole in eq. (5.24) with a specific diagram, as we have done in figs. 4 and 8. 
Although the definition of the completed-cut terms C„ is not unique, the ambiguity cancels 
between (^^(O), Inf C^, and the sum over CRn residues in 0„. 

In general, it is useful to combine all the rational functions not included in the cut 
completion into a single function, 

Rn = — Inf An - Inf C„ + i?^ + 0„ . (5.30) 

Cr 

so that 

An(0)=cr[c„(0) + P„] . (5.31) 

If C„(0) is chosen to preserve a symmetry of the amplitude {e.g., under a particular permu- 
tation of legs), then i?„ will also have this symmetry, even if the individual components of 
Rn do not. 

D. Determining Terms Arising from Large Shifts 

Our remaining task is to determine the large-z contributions set aside in eq. (5.5), and 
put back unevaluated in eq. (5.29). As already discussed in sections III and IV, to do so we 
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FIG. 10: A non-standard factorization using an auxiliary [a, h) shift is assumed to be suppressed 
in the large-z hniit of the primary [j, I) shift, if the tree amphtude naively appearing in the factor- 
ization contains both legs j and / and is suppressed. 

will use a second, auxiliary shift, 

[a,6) : Aa^Aa-wAfe, Xh^Xh + wXa, (5.32) 

distinct from the primary shift in eq. (5.1). It is useful to choose the auxiliary shift so that 
Aniw) vanishes in the large w limit. (For the gluon amplitudes we consider in this paper, 
and likely in general, such choices can be found, as we discuss in section VI.) The price 
that we must pay is the presence of contributions to the amplitude from channels with non- 
standard complex factorizations. However, we can arrange matters so that these channels 
do not contribute to the terms we are seeking to compute, those that survive in the limit 
when the original shift's parameter z becomes large. 

To do that, it suffices to ensure that channels with non-standard complex factorizations 
vanish at large z when the primary shift (5.1) is applied to the auxiliary recursion relation. 
The details of the factorization behavior in those channels are then unimportant. 

For multi-particle channels, factorization in complex momenta is the same as in real mo- 
menta. For two-particle channels with opposite-helicity gluons, the discussion in section III B 
showed that only the standard tree-level tensor structure contributes, and that for scalars in 
the loop the relevant form factor actually vanishes. For like-helicity gluons, in the complex 
kinematics for which the tree vertex is non- vanishing, the situation is more complicated. The 
structure of the factorization is known empirically for one of the two possible helicities of 
the intermediate gluon, for the case where the other possible intermediate helicity vanishes. 
This case occurs in studying the finite series of n-gluon amplitudes A^~^°°p(1~, 2+, . . . , n"*"). 
The factorization has the form [30], 



A^-i-P(l,...,6+,(6+l) 



,n) 



vl-'^^^{-K\ b\ {b + ir)A':^,{l, ...,K-,...,n) (5.33) 
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X (1 + correction factors) . 

This configuration is shown schematically in fig. 10(a). So long as the legs shifted under the 
original, primary, shift are both contained on the A^'^^'^ side of the factorization, the form of 
the correction factors in eq. (5.33) is such that the term vanishes as the primary shift variable 
z -^ oo, because the tree amplitude vanishes in that limit. We will defer a general study 
of the structure for the other intermediate helicity shown in fig. 10(b); the only property 
that we will need is the vanishing of the primary shift z ^ oo limit for contributions where 
both shifted legs are on the tree side of the factorization, as in fig. 10. This property holds 
for previously-computed five- and six-point amplitudes [27, 48] , and we will assume it holds 
more generally for both intermediate helicity configurations. The validity of this assumption 
can be tested at the end of a calculation, by checking all the symmetry and factorization 
properties of the fully assembled amplitude. 

The auxiliary shift gives us the following expression for the complete amplitude, 

An{0) = cr \Cn{0) - Inf dn + i?^'["''> + Ol:M , (5.34) 

L [a,b) J 

where the superscript [a, b) denotes the legs shifted under the auxiliary shift, and where the 
recursive diagrams are built and the overlap contributions determined, with respect to this 
shift. The large-parameter (large-w) terms of the auxiliary shift (5.32) are absent by design. 
We can now extract the large-parameter (large-z) behavior with respect to the primary [j, I) 
shift, 



Inf An = Cr Inf 



d.(0)-InfC„ + i?^'I"'^> + O^^> . (5.35) 

[a,b) J 



Following the discussion above, we arrange the shifts so only channels with standard factor- 
izations will survive in the laige-z limit of the [j, I) shift. In extracting the laige-z behavior, 
we must in general keep all non-vanishing contributions to the amplitude, which may arise 
not only from the recursive diagrams with respect to the auxiliary shift, but also from 
completed-cut or overlap terms. In many practical cases, however, the only surviving con- 
tributions are from a limited set of recursive diagrams. Indeed, the typical surviving term 
will have the form, 

lniR{kp„...,ka,...,kp_,,-P'')x^xA'''''%kp^,...,h,...,kp ^,P-''), (5.36) 

where legs j and / are both on the loop side. 
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FIG. 11: For a [j,l) shift, the four potential channels with non-standard complex singularities. 
Whether these are actually present in a given amplitude depends on the helicities of the legs nearest 
to j and / in the color ordering. 

In the next section, we shall present shift choices for general helicity configurations that 
implement the approach outlined in this section: a primary shift free of non-standard chan- 
nels, but having non-trivial large shift-parameter behavior, and an auxiliary shift free of 
non-trivial large shift-parameter behavior, but containing non-standard channels that in 
turn vanish at large values of the primary parameter. 

VI. GENERAL HELICITIES 

As an illustration of our strategy, in this section we now present specific shift choices 
for determining the rational-function parts of generic one-loop n-gluon amplitudes. As 
discussed in previous sections we must choose a primary shift so that non-standard complex 
singularities do not occur in the recursion. If there are contributions from large values of 
the shift parameter z, we determine these using an auxiliary shift and recursion relation. 

A. Empirical Structure of the Amplitudes 

In order to proceed, we will need a few analytic properties of amplitudes. Unfortunately, 
as yet there are no theorems to guide us on the properties of complex factorizations of 
amplitudes or on the behavior of loop amplitudes under large complex shifts z. We therefore 
follow the empirical approach of refs. [27, 30, 31]. We observe certain useful properties of 
known amplitudes and then use these properties to aid in the computation of new amplitudes. 
We shall not prove these properties, noting that such proofs would be valuable to help guide 
future developments. This empirical approach has been effective for obtaining a variety of 
new one-loop amplitudes [27, 30, 31, 32]. By now, a large number of QCD amplitudes are 
known analytically [30, 31, 32, 47, 48, 57, 63, 64, 65, 66, 67], making it straightforward to 
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develop a heuristic understanding of their analytic properties. 

Our confidence in this pragmatic approach stems from the rather non-trivial checks that 
may be performed on any loop amplitude. In particular, the factorization properties of 
one-loop amplitudes in real momenta are well understood [5, 6, 59, 60, 61] and provide 
rather non-trivial constraints on the amplitudes by demanding that every pole in the ampli- 
tude corresponds to a physical factorization. The non-trivial consistency requirement stems 
from the fact that only a limited number of factorization channels enter into the recursive 
construction. 

An investigation of the analytic properties of the known one-loop amplitudes reveals some 
striking properties: 

1. For any [—, +) shift, a//n-gluon amplitudes vanish for a large shift parameter, z -^ cxd. 

2. For n-gluon amplitudes with the split helicity configuration, A„(l~, 2~, . . . , (m — 
1)^, -m^, (m + 1)^, . . . , n"^), an alternative set of shifts where the amplitudes vanish at 
large z are [l,m — 1), [m — 1, 1), [m,n) and [n,m). 

3. For a given [j, I) shift, there are no more than four channels with non-standard complex 
singularities, depending on the helicities of the legs nearest to j and /, as depicted in 
fig. 11. 

The above properties are by no means exhaustive. In particular, there are other shifts 
where one-loop n-gluon amplitudes vanish at large z, although the above observations will 
be sufficient for our purposes here. 

As already discussed in section VD, an empirical rule for suppressing diagrams with 
non-standard complex singularities in an auxiliary recursion relation is to ensure that the 
primary shift legs are both on the tree side of the naive factorization and that this tree 
amplitude is suppressed in the large-z limit of the primary shift. Such configurations are 
displayed in fig. 10. 

It is worth mentioning that based on our empirical studies of MHV supersymmetric 
amplitudes [5, 6, 11], it appears that in the supersymmetric case, the complete set of shifts 
where An{z) vanishes for 2; ^ oo is identical to the set of shifts where this is true at tree 
level, i.e. any [—,+),[+,+) and [— , — ) shift. 
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B. Systematics for General Helicities 

Using the above empirical observations, we now present a systematic procedure for finding 
pairs of shifts which will allow us to compute the rational terms in any n-gluon amplitude 
while avoiding non-standard complex singularities. 

Depending on the helicity configuration, we will use the three independent shift choices: 

• If the amplitude contains four color-neighboring legs having the helicity structure 
i+, {i + 1)-, {i + 2)+, {i + 3)-, then choose the shift [i + 1, i + 2). With this shift, the 
amplitude vanishes as 2; ^ oo and no non-standard complex singularities appear in 
the recursion relation. Only a single shift, and hence only a single recursion relation, 
is required in this case. 

• If the amplitude has three nearest neighboring legs i~, (i + l)~, (i + 2)"'", choose [i + 1, i) 
as the primary shift. For determining the behavior of the amplitudes for large values of 
the primary shift parameter, choose an auxiliary shift [a, h) such that a is a negative- 
helicity leg, 6 is a positive-helicity leg and a ^ i — l^i^i + 1. 

• For the special case of split helicity configurations, A„(l~, 2~, . . . , (m — 1)~, -m^, {m + 
l)"*", . . . , n"''), a rather convenient choice is a primary [1, 2) shift and an auxiliary [n, m) 
shift. 

The above choices are not the complete set of choices that we need. However, all the 
remaining cases are simply related to the above ones via parity conjugation or a reversal of 
legs in the color ordering. For convenience, we also list these shift choices: 

• If four neighboring legs in the color ordering have helicities i~, (i + 1)"^, (i + 2)~, (i + 3)"'" 
then choose the single shift [i + 2, i + 1). 

• If the amplitude has three nearest-neighboring legs «'*',(« + 1)~, (^ + 2)^, choose 
[i + 1, i + 2) as the primary shift. As the auxiliary shift choose any [a, h) such that a 
is a negative- helicity leg, 6 is a positive-helicity leg and a7^i + l,i + 2,i-|-3. 

• If the amplitude has three nearest-neighboring legs i~, {i + 1)"^, {i + 2)"*" choose 
[i + 2, i + 1) as the primary shift. As the auxiliary shift choose any [a, h) such that a 
is a negative-helicity leg, 6 is a positive-helicity leg and 67^i + l,i + 2,-i + 3. 
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TABLE I: This table lists helicities and pairs of shifts that may be used to construct any six- or 
seven- gluon amplitude with three negative- helicity legs. The primary shift in the second column 
generates a recursion relation which does not have non-standard complex singularities, but the 
amplitude may not vanish for a large shift parameter z. Under the auxiliary shift in the third 
column, the recursion relation may have non-standard complex singularities in the channels listed 
in the fourth column. However, these singularities should be suppressed in the large-z limit of the 
primary shift. 



Helicity 


Primary shift 


Auxiliary shift 


Suppressed channels 


l-,2-,3-,4+,5+ 


[1,2) 


[5,4) 


— 


l-,2-,3+,4-,5+ 


[4,3) 


— 


— 


l-,2-,3-,4+,5+,6+ 


[1,2) 


[6,4) 


(45) 


l-,2-,3+,4-,5+,6+ 


[4,3) 


— 


— 


l-,2+,3-,4+,5-,6+ 


[1,2) 


— 


— 


l-,2-,3-,4+,5+,6+,7+ 


[1,2) 


[7,4) 


(45) 


l-,2-,3+,4-,5+,6+,7+ 


[4,3) 


— 


— 


l-,2-,3+,4+,5-,6+,7+ 


[1,2) 


[5,6) 


(67) 


l-,2+,3-,4+,5-,6+,7+ 


[3,4) 


— 


— 



With these choices we should then be able to construct the rational-function contributions 
of all unknown n-gluon amplitudes, once the cut-containing pieces are known. (The above 
choices are not useful for constructing amplitudes with identical helicities, but those are 
already known [66, 67].) If more than one of the above choices is satisfied in a given 
amplitude, one may choose whichever is the most convenient. In Table I we have listed 
all the helicity configurations with three negative-helicity legs and up to seven external 
gluons, along with choices of primary and auxiliary shifts which may be used to construct 
the amplitudes. (In the first row, for 1^, 2^, 3^, 4+, 5+, our choice of auxiliary shift actually 
has no non-standard factorization channels, so it could be used by itself to fully determine 
the amplitude. We display this particular shift choice because it is based on the above rules 
and generalizes to the case of more adjacent positive-helicity gluons.) 

It is important to note that there are many other valid shift pairs besides those in the 
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TABLE II: This table lists examples of valid shift pairs besides those selected by the procedure 
discussed in the text. 



Helicity 


Primary shift 


Auxiliary shift 


Suppressed channels 


l+,2+,...,n+ 


[5,7) 


[1,3) 


— 


l-,2+,...,i-,(j + l)+,...,n+ 


[1,J) 


— 


— 


l-,2-,3-,4+,...,n+ 


[1,2) 


[4,5) 


(56) 


l",...,(m- l)",m+,...,n+ 


[1,2) 


[m,m + 1) 


(m+l,m + 2) 



above construction. For example, although we can use the rules to determine all amplitudes 
with two negative helicities, it turns out that a somewhat more convenient choice is to 
choose a [i,j) shift, where legs i and j are the two negative-helicity legs, as we shall discuss 
in a companion paper [47]. In many cases, it is also possible to relax the conditions we 
impose on the shifts. For example, we have been demanding that under the auxiliary shift 
the amplitude vanish for large shift parameter. In fact, this restriction is not necessary; 
we need only demand that any such terms do not contribute to the large shift terms of 
the primary shift. For example, the identical-helicity amplitudes A„.i(l+, 2+, . . . ,n''') can 
be determined using a primary [5,7) shift and an auxiliary [1,3) shift (for n > 7), even 
though the amplitudes do not vanish [30] under the large- 2; limit of either shift. In Table II 
we have collected a variety of examples of shift pairs which may be used to determine 
the amplitudes, but are outside the class of shifts described above for determining general 
helicity configurations. 

We expect that a similar strategy will be effective for amplitudes with massless quarks, 
and for amplitudes with external massive vector bosons or Higgs particles. With suitable 
modifications it should be possible to use the on-shell bootstrap to construct amplitudes 
with massive particles in the loops as well. 

VII. RECURSIVE DETERMINATIONS OF LARGE-z BEHAVIOR 

Following the procedure of the last section we now determine the large shift-parameter 
behavior of some sample amplitudes. We focus on amplitudes with three or four color- 
adjacent negative-helicity legs, as a non-trivial illustration of the method. Because the 
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logarithmic terms in these amphtudes have already been calculated [26], we can obtain the 
complete amplitudes by computing the rational terms recursively, as we do in the next 
section. We can confirm indirectly that our approach to determining the large-z behavior is 
valid in these cases, by verifying that the amplitudes have the proper symmetries, and that 
they factorize correctly. 

To determine the large-z behavior of an amplitude, with the help of an auxiliary shift, 
one must include in general the large-z behavior of the completed-cut, overlap and recursive 
contributions (with respect to the auxiliary shift). In subsection C, we shall encounter a five- 
point example where all three types of contributions are nonvanishing at large z. However, 
in subsections A and B, we shall arrange the shifts, for the respective cases of three and 
four color- adjacent negative helicities, so that the entire large- 2; behavior comes from the 
recursive diagrams of the auxiliary shift. In this way it is very simple to obtain compact 
expressions for the large- 2; behavior of these amplitudes under the [1, 2) shift. 



A. <r0(l-,2-,3-,4+,5+,..., 



n ' 



Consider now the n-gluon amplitudes with three color- adjacent negative-helicity legs. We 
use the special split-helicity shift choice described in the last section: a [1,2) shift as our 
primary shift, and an [n, 4) shift as the auxiliary shift. This shift allows for an especially 
simple determination of the large-z behavior in the n-point case. 

How do the individual contributions behave as the z parameter associated with a [1,2) 
shift is taken to infinity? First consider the completed-cut contributions, C„. These contri- 
butions were obtained in ref. [26] and are given in section VIII, eq. (8.1), where we describe 
them in more detail. As explained in that section, Cn is non-vanishing in the large-z limit 
of the [1,2) shift (see eq. (8.10)). However, a computation shows that the large-z limit of 
the overlap terms obtained from C„ cancels this term. So there are no laige-z [1, 2) shift 
contributions arising from the sum of the completed-cut and overlap contributions of our 
auxiliary [n, 4) shift. All that remains is to inspect the large-z behavior of the recursive 
diagrams. 

As a warm-up for the all-n case, we consider the six-point case again, this time with an 
auxiliary [6,4) shift instead of [3,4). The recursive diagrams are shown in fig. 12, omitting 
diagrams where the tree vertex vanishes. The simplest diagrams to analyze are those where 
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FIG. 12: Recursive diagrams for ^^f °(1", 2", 3",4+, 5+,6+), for the auxiliary [6,4) shift. We 
drop the diagrams for which the tree vertex vanishes. Diagram (e) has non-standard complex 
singularities. Only diagrams (k) and (1) contribute to the large-z limit of the [1,2) shift. Diagrams 
(c) and (f) represent the factorization-function contribution. 

both shifted legs 1 and 2 are attached to a tree vertex in a 'standard' channel, as is true 
for figs. 12(a)-(c). The remainder of the diagram (the propagator and other vertex) is 
independent of the shift parameter, and so the large-parameter behavior is determined by 
the tree alone. Because this shift is a well-behaved shift at tree level [29, 44], the tree 
vertex, and hence the diagram, vanishes in the large-z limit. The same reasoning applies 
to diagrams where both shifted legs are attached to a loop vertex, and where we already 
know from previous computations that the [1, 2) shift is well-behaved. This is the case 
for fig. 12(d). We cannot be certain that the same logic will apply to diagrams such as 
fig. 12(e), where both shifted legs are attached to a tree vertex, when the contribution is in 
a channel with non-standard factorization. Nonetheless, analysis of specific analogous cases 
leads us to conclude that these diagrams should be suppressed in the large-parameter limit. 
As discussed in the previous section, this assumption can be tested once the final answer for 
the amplitude is in hand. Diagrams where the shifted legs are attached to different vertices, 
such as those in figs. 12(f)-(j), must be analyzed explicitly. These diagrams involve only 
standard channels, and all do indeed vanish in the large-parameter limit. 

This leaves us with two diagrams, those of figs. 12(k) and (1). Here, the shifted legs 
are both attached to the same five-point loop vertex with a standard factorization, so the 



47 



large-parameter behavior is determined solely by the loop vertex. We know from eq. (3.14) 
that it has non-trivial behavior in that limit, which survives to contribute to the auxiliary 
recursion relation. Collecting these two contributions, we obtain the desired relation. 



\Af=0/T- o- Q- 



InfA^r(l-,2-,3-,4+,5+,6- 



Inf <r°(l-, 2-, ku, 5+, 6+) — ^r (3-, 4+, -Kt,) (7.1) 

[1>2> S34 

+ Inf A^=\r, 2-, 3-, i^+, 6+) — ^r (4+, 5+, -K^,) , 

11.2) S45 

where a refers to legs shifted and frozen according to the auxiliary [6,4) shift. Evaluating 
this relation (and as we shall see below, even solving it for all n) is straightforward, thanks 
to its similarity to a tree- level MHV recursion relation [29]. Plugging in the known values 
on the right-hand side, we obtain, 

.^=on- o- .-.+ .+ «+._ ^^r (li^34)'(26> 1 [k,Af 



Inf A^f°(r,2-,3~,4+,5+,6^ ., „ . . „ 

[1,2) 6,1 V , , , , , . 3 ^16^2 ^^^^^^ ^5g^j^2] S34 [34] [7^34 3] 

icr (13)^(2 6) 1 [5 4]' 



+ 



3 (16)' (3i^45) (i^45 6)[12] S45 [7^45 5] [47^45] 
^ tcr (13)^(2 6) 

3 (16)' (3 4) (4 5) (5 6) [12]' ^ ' ' 

in agreement with eq. (4.4). 

It is straightforward to generalize eq. (7.1) to obtain the large-z behavior of the n-gluon 

amplitude under a [1,2) shift. Following the same logic, we obtain the recursion relation, 

using an auxiliary [n, 4) shift, 

Inf<r°(l",2-,3-,4+,5+,...,n+) = 

[li2) 

Inf At=Ul-, 2-, i^3-4, 5+, . . . , n+) — ^r (3-, 4+, -K^,) (7.3) 

11.2) S34 

+ Inf Alti(l", 2-, 3-, i^+ , 6+, . . . , n+) — Af^ (4+, 5+, -K,,) , 

|1.2) S45 

where we assumed that diagrams involving non-standard factorization, and diagrams where 
legs 1 and 2 appear on different vertices, are suppressed at large z, the same way they were 
at six points. Solving this recursion relation by induction, we obtain the large z behavior at 
n points, 

iAr=o/i- o- Q- /1+ ^+^_^cr (13)^(2n) 



lniA^=\l-,2~,3-,4\.--,n^) = -^ ^ Z I ' • (7-4) 

[12](nl)'n(M^ + l)) 



[1,2) "'' ^ ' ' ' ' ' ' 3 

[izj^ni; 

A;=3 
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FIG. 13: Non-vanishing recursive diagrams for A^^^[l^ , 2^, 3^, 4^, 5+, 6+), for the auxihary [4, 5) 
shift. We use this recursion to determine the large-z behavior under the primary [1, 2) shift. Only 
diagrams (a) and (b) contribute to the large-z hmit. 

We will use this result in the next section to construct a recursive solution for the rational 
functions in the n-point amplitude. 

B. ^^f0(l-,2-,3-,4-,5+,...,n+). 

We now consider the case of four color-adjacent negative helicities. Once again it is useful 
to illustrate the six-point case, Agf='^(l~, 2~, 3~,4~, 5+, 6+), before turning to the n-point 
case. Following the discussion of the previous section, a convenient pair of shifts is the 
[1, 2) shift as the primary shift and the [4, 5) shift as the auxiliary shift for determining the 
large-z behavior of the amplitude under the primary shift. (An alternative choice would be 
an auxiliary [n, 5) shift.) 

As was the case for three color-adjacent negative helicities, an examination of the cut 
terms given in ref. [26] reveals that in the large-z limit of the [1, 2) shift there are no 
contributions from the sum of the completed-cut and overlap contributions. So here we 
focus on the non-vanishing rational-recursive contributions. 

The auxiliary [4, 5) shift diagrams are shown in fig. 13 and are given by 

i?f''''>(l-,2-,3-,4-,5+,6+) = /)(-) + D(b) + Z}W + D(d) + /)(e) 

= Af'{3', 4-, -K+,)—R,{l', 2-, K,„ 5+, 6+) 

+ R,{1', 2-, 3-, 4-, i^5+)— ^r (5+, 6+, -K^,) 

S56 
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+ i?3(3-, 4-, -kt,)—Af'^{r, 2-, k,,, 5+, 6+) 

S34 

+ i?3(3-, 4-, -k^^)—Af\l~, 2-, i^3+ , 5+, 6+) 
S34 

+ i?4(2-, 3-, 4-, -i^+J^^r (1-, i^234, 5+, 6+) , (7.5) 

S234 

where a denotes momenta shifted and frozen according to the auxihary [4, 5) shift. We apply 

the [1, 2) shift to this recursion relation and then take the large- 2; limit of this primary shift. 

We assume that diagrams (c) and (d), which contain non-standard complex singularities, 

are suppressed in this limit (as discussed in section VD), because the corresponding tree 

vertices are suppressed. 

Diagram (e) is also suppressed, but that requires a closer inspection. To analyze diagram 

(e) we need the vertex, obtained by parity conjugation from eq. (A14), 

o n+ 2- V 4-^ - ' [2 4] (2 4)^ 

^'^^ '2 '^ '^ )-3(12)[23][34](41)- ^^'^^ 

Then diagram (e) is. 



D(^) = i?4(2-,3-,4-,-i^+34) Ar(l-,i^234,5+,6^ 

■§234 



,3 



3 



t [24] (2 4) 1 (li^234) 

3 (i^234 2)[2 3] [34] (4i^234)S234 (ii'234 5) (5 6) (61) 

I [24] (24)^ 1 (1-|^234|5-) 



-\3 



(7.7) 



3(2-|^234|5-)[2 3] [3 4] (4- 1 ^234 1 5- )S234 (5- 1 ^234 1 5-) (5 6) (61) 
We can simplify this term further, by substituting in the remaining hatted variables, which 
become independent of 2; as 2; ^ oo; but already in this form we can see that it vanishes 
in the large-z limit of the [1, 2) shift. Therefore diagram (e) does not contribute to the 
recursion relation for the large-z terms. 

Thus we obtain a very simple recursion relation for the non-vanishing large-z terms, based 
only on diagrams (a) and (b), 

Inf <r°(r, 2-, 3-, 4-, 5+, 6+) = Ar(3-,4-,-i^3+)— Inf<r°(l-,2-,i:34,5+,6+) 

[1.2) S34 [1,2) 

+ ^r (5+, 6+, -k^,) — Inf <r°(l", 2-, 3-, 4-, i^+) • 

(7.8) 

On the right-hand side there are two large-2; contributions. The first one is familiar from 
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eq. (3.24), 



Inf A^f°(l-,2-,3-,4+,5- 
[l,2> 5,1 V ' 



ct 



(13)' (2 5) 



(7.9) 



3 [12] (5 1)^(3 4) (4 5)- 

To get the second one, we start from formula (A16) for A^j"^ (1^,2^, 3^,4'^, 5^), relabel, 
and take the parity conjugate, to get. 



<r°(l-,2-,3-,4-,5- 



cr 1 



(14)= 



[5 3]^ (3 4) [2 4] [5 2]^ (21) [31] 



3 [2 3]' L (51) (4 5) [51] [12] [3 4]' [5 4] [4 3] [2 1]' 

(7.10) 

Applying the [1,2) shift and extracting the z^ terms gives us the desired contribution, 

which we wish to feed into the recursion. 



Inf A^f°(l~,2-,3~,4-,5" 

[1,2> 5,1 V ) 



cr 1 



(14)' [5 2]' (21) [31] 

(51) (4 5) [5 4] [4 3] [21]' 



(7.11) 



3 [2 3]' 

Relabeling and inserting the large- 2; terms (7.9) and (7.11) into the recursion (7.8) we 
obtain, 

Inf<r°(l-,2-,3-,4-,5+,6+) 

= ^r (3-, 4-, -i^3+) — Inf <r°(l", 2-, k,„ 5+, 6+) 

S34 11,2) 

+ Ar%5\ 6+, -k^,) — Inf <r°(i-, 2-, 3-, 4", i:+ ) 

S56 [1:2) 



Cr 



(34) = 



(1^:34)' (2 6) 



3 (47^34) (K343) S34 [12] (61)' (i^34 5)(5 6) 



+ i 



Cr 



[5 6]= 



1 1 



3 [5 7^56] [^56 6] S56 [2 3]^ 



y^y 



[i^56 2]'(2 1)[31] 

(i^sel) (4i^56) [i^56 4] [4 3] [21]' 



Cr 



-\3 



(r|(3 + 4)|5-)''(2 6) 



[12][3 4][4 5](61)'s345(6-|(4 + 5)|3-) 



(7.12) 



+ 



(14)= 



^ (4^1(5 + 6)12')' (21) [13] 

[2 3]' (4 5) (5 6) (6 1) [1 2]' [2 3]' (4 5) (5 6) S456 (6-| (4 + 5) [3") 



We have checked that this result is in agreement with the large-z behavior extracted from 
the parity conjugation of the six-point amplitude obtained in ref. [27]. 

It is straightforward to generalize the discussion to n-point amplitudes. Following the 
same logic as for six points, we obtain a recursion relation for the large- 2; behavior under 
the [1,2) shift. 



Inf<r°(l",2-,3-,4-,5+, 

|1,2) 



,n 



(7.13) 
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FIG. 14: Nonvanishing diagrams of the [4, 5) shift recursion relation for A^-^ ^{1 ,2 ,3 ,4+, 5'^). 
Only diagram (a), which has legs 1 and 2 on opposite sides of the pole, contributes to the large-z 
limit of the [1,2) shift. 



^ ^ " .Af=o 



^r (3-, 4-, -7^3+ )— Inf Alti(l-, 2-, K~„ b\...,n' 

■S34 11.2) 



+ Ar%h\ Q\ -K^e) — Inf ^lti(l"' 2", 3", 4", K^„ 7+ 



, . . . , 



n 



We have solved this recursion relation with the result, 



InfA^r°(l",2-,3-,4-,5+,...,n+) (7.14) 

|J-i2) 



.cr 
t — 
3 



(14)' , (4-|^5...j2-)='(21)[13] 



+ 



n— 1 n— 1 

L [2 3]' n (^ (^ + 1)) (^ 1) S4...„ (n- 1 ^4...„ |3-) [1 2]' [2 3]' U{k{k + 1)) 

fc=4 fc=4 

, Y^ (j(j + l))(l'|^3....^5...,|4+)'(2n) ■ 

"*" 2-^ n-l 

^=5 S3...JS,.., iu + 1)-| ^4..., |3-) {r\^4...j |3-) [1 2] n (^ (^ + 1)) (n if 

k=A 

C. Another Look at A^f °(1-, 2-,3",4+, 5+) 

In the two previous subsections, the completed-cut and overlap contributions did not 
contribute to the large-z terms of the primary shift. However, as we now illustrate, 
this is not always true. We reexamine the large-z limit of the five-point amplitude 
A^j"'^(l^, 2^, 3^, 4+, 5^) under a [1,2) shift. This example was already considered in sec- 
tion III. However, here we recursively construct the large-z behavior by using a different 
auxiliary shift, a [4, 5) shift. Although there are no non-standard factorization channels here, 
the example displays several important features: The large-z recursive contribution comes 
entirely from a diagram where legs 1 and 2, the primary shift legs, are split across the pole 
of the auxiliary shift. Also, as was just mentioned, we need to account for contributions 
from completed-cut and overlap terms under the auxiliary shift. 

The nonvanishing recursive diagrams of the auxiliary [4, 5) shift are shown in fig. 14. 
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First observe that the behavior of diagram (b) under the [1,2) shift is determined by the 
rational part of the four-point one-loop MHV amplitude with adjacent negative helicities, 
which is proportional to the tree amplitude, according to eq. (A21). Tree amplitudes vanish 
at large z for shifts of identical-helicity pairs of legs [29, 44]; hence diagram (b) does not 
contribute to the large-z limit. Therefore, we need only evaluate the contribution of diagram 
(a). Here we find a new feature: The shifted legs of the primary [1, 2) shift cross the pole 
in the auxiliary recursion relation. Nevertheless, the large- 2; contribution can be determined 
just as in the previous examples. This pole-crossing contribution is given by, 

S51 
_ i [i^5i3] (i:5i3)' 1 [5i^5l]' .7^3. 

3 [2 3] (3 4) (4 i^5i) [i^5i 2] S51 [1 i^5i] [51]' 
where a refers to legs shifted according to the auxiliary [4, 5) shift. We then obtain in the 
large-^ limit of the [1,2) shift the contribution, 

Infi?f'l^'^>(l-,2-,3-,4+,5+) = -- ^^^M . (7.16) 

[i,2> '' V , , , , ; 3 [2 3] (3 4)' (4 5) (51) ^ ' 

Comparing eq. (7.16) to the known expression in eq. (3.24) we see that the two expressions 

do not match. However, we also have to take into account the contributions to the large-^ 

limit from the overlap contributions as indicated in eq. (5.29). We apply the [4,5) shift to 

Ci?5, which can be extracted from eq. (3.16). We obtain the overlap contribution from the 

residues of CR^{w)/w at the following values of w, illustrated in fig. 15, 

{14)' |35] ^ ' 

Evaluating these in the standard way, we obtain the large-z contributions, 

I^f 0(a), [4,5) ^ _^^13)^]25]j23^ _ I (13)^ [4 5] 

[i,2> ' 3 [12] [2 3] (3 4)' (51) 6 [1 2] [23] (34) (5 1) ' ^' ^ 

InfOf)'[^'^> = liipi. (7.19) 

[i,2> ' 6[12][2 3]S34 ^ ' 

Notice that both eqs. (7.16) and (7.18) have double poles in 1/ (3 4), which cancel correctly 
upon adding up all contributions. 

Finally, we need to take into account the contribution from the completed-cut terms C5, 
which are easily found to be, 

Infa.(l-,2-,3-,4^5-)^-lii^IMilMl^m4^(iM^. (r.20) 

[1,2) 5V , , , , ; g [12][2 3](51)'s34 
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FIG. 15: The overlap diagrams of the [4, 5) shift. 



Upon adding all contributions, eqs. (7.16), (7.18), (7.19), and (7.20), we obtain 



j\f=o 



Inf A^.l 

[1,2> ^'' 



-ICy 



1 



(13)^(2 4) 



1 (13)^ [2 5] (2 3) 



(13)^ [4 5] 



3 [2 3] (3 4)' (4 5) (51) 3 [1 2] [2 3] (3 4)' (5 1) 6 [1 2] [2 3] (3 4) (5 1) 
1 (13) [4 5]^ 



cr 



6 [1 2] [2 3] S34 6 
(13)^^(2 5) 



^ 1(1 2) [2 4] (13) ((12) [2 4] -(13) [3 4]) 



[12][2 3](51)"s34 



(7.21) 



3 (15)' (3 4) (4 5) [12]' 

in agreement with eq. (3.24). 

This example demonstrates that, in general, if we wish to obtain the correct large-^ 
behavior using an auxiliary recursion, we must include the contributions from the completed- 
cut and overlap terms. In practice, however, often only the recursive diagrams contribute. 



VIII. COMPLETE SPLIT HELICITY AMPLITUDES 

In section VII A we determined the large-z behavior, under a [1,2) shift, of the n-point 
amplitudes with three nearest-neighboring negative helicities, A^f °(1~, 2~, 3~, 4^, . . . , n^). 
We now use this result to evaluate the remaining rational terms, using the formalism de- 
scribed in section V. The n-point tree amplitudes for these helicity configurations were 
determined in refs. [9, 34, 43, 68], and are given in eq. (A9). 

A. Cut Contributions 



The completed-cut contributions for general split-helicity A/" = amplitudes have been 
obtained in ref. [26]. For the case with three negative helicities, the cut part of the amplitude 

<r°(l-,2-,3-,4+,...,n+)is, 

1 



C„(l-,2-,3-,4+,..., 



n'j = -^<r'(l",2-,3-,4+,...,n^ 
3cr 
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+ -Ar(l-,2-,3-,4+,...,n^ 



n—l 
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o / J "'n,r 



L2(^^^ 



W=0 ^V-S2 



„3 

^2...r 



n-2 
r=A 



u<^^ 



"^ ■ '^2...(r+l) 



r=4 



'^S-.-Cr + l) ■ 



'3...(r+l) 



l.l) 



where 



d: 



9n 



W=0 ^ (3-| ^3...,4 |1 + ) (31 fe^2...r |1 + ) (3-| ^3...rfe, ^2...r]^2...r |1 + ) 

r— 1 n 

k=3 k=r+l 

M=0 _ ^ (3-|^3...0-+3)^2...0-+3)[^r+l,^2.J|l + )((j+3)(j+4)) 



f8.2) 
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X 

S3...{i+3) S2...(i+3) 
>< (3 I -^3...(i+3)-^2...(i+3)-^(r+l)...ll^r+l |l ) 



h 



M=0 



^J yn,n-r+2 



(123...n)-^{321ra...4) 

Here ^^f ^ is the contribution of an A/" = 1 supersymmetric chiral multiplet consisting 
scalar and a fermion running in the loop. The result for this amplitude is [26, 49], 

<fi(l-,2-,3-,4+,...,n+) 



(8.3) 
(8.4) 

of a 



A, 



tree 



cr -^ (Ko(s„i) + Ko(s34)) -T^crJ^dn, 
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2' 



n-2 Ln( r'"-'- ) .• !^2 Lo( ~'^ ■'• 



9n,r 
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■52...(r+l) 2 
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S2...r 
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■53...{r + l) ■ 



■53...(r+l) 
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where, 
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r— 1 n 
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FIG. 16: The recursive diagrams of the [1,2) shift for an n-point hehcity amphtude. Some of the 
diagrams may vanish depending on the hehcity choices. 

From eq. (8.1) we can extract the rational parts of the cut completion, Ci?„. These terms 
are given by, 

CRn = (^^ + 0Ar(l-,2-,3-,4+,...,n+) 

. n— 1 . n— 2 

_ ^ V^ ,7-^=0 '^3...r + S2...r ?_ V^ -A/'=0 '^2...r + ■52...{r+l) 



. S2...rS3...r{S3...r ~ S2...r) . S2...rS2...{r+l){S2...(r+l) — S2...r) 



. n-2 

- ^E'e° '"V'""" ^' (8^9) 

6^ ■^3...rS3...{r+l)(S3...(r+l) - S3...r)^ 

This cut completion has a (constant) boundary contribution as z ^ oo under a [1, 2) shift, 
stemming from the d'-'^'^^^-term. It is given by, 

,r^ ^^ (3-|^3....fe|l+) (13) (g3...r(13) + (3l^3....^2|l+))(r(r + l)) ,^ ^^, 

''^ -4 S3...r(2+|^2...r|r+)(2+|^2...rl(r + in n(MA; + l)) 

A:=3 

B. Recursive Contributions 

We now discuss the recursive diagrams of the [1,2) shift, shown in fig. 16, for a generic 
helicity configuration. As explained in section V, we must sum over both helicities of the 
gluon crossing the pole, and over loop versus tree vertices appearing on either side of the pole, 
as well as factorization-function contributions [59] in multi-particle channels. For specific 
helicity configurations some of the diagrams may vanish. 

For the case of A^f °(1~, 2~, 3~, 4"*", . . . , n"*"), the recursive contributions are. 



C(l-,2-,3-,4+,...,n^ 
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(r+l)+ 




FIG. 17: Overlap diagrams corresponding to the S2...j, -52. ..r and S2...(r+i) channels. 
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j=4 



S2.. 



xAf-(2-,3-,4+,...,^+,-i^. 
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2...i) 



Y,Rn-^+2{l', Kt,, {i + 1)+, . . . ,n+)^Ar^(2-,3-,4+, . . . ,^+, -i:. 



2...i> 



fS.lll 



where the last Une comes from the factorization-function contributions. We refrain from 
quoting the exphcit expression for R^ here. It can be obtained straightforwardly by inserting 
the known tree amplitudes [33, 43, 69] (quoted in the appendix) and rational parts of loop 
amplitudes [31, 32] into eq. (8.11). 

An interesting feature of this recursion relation is that all amplitudes on the right-hand 
side have fewer than three negative helicities. In contrast, the recursion relation found for 
two negative-helicity gluons [27] contains on the right-hand side lower-point loop amplitudes 
with the same number of negative helicities as on the left-hand side, namely two. For this 
reason, to solve that recursion relation in closed form for all n required an 'unwinding' 
procedure [32]. Unwinding is not necessary in the present case because all terms on the 
right-hand side are known — given the two- negative-helicity solution of ref. [32]. 

C. Overlap Contribution 

Finally, we need to compute the overlap contribution to avoid double counting between 
the recursive diagrams, eq. (8.11), and the completed-cut terms, eq. (8.9). First we observe 
that the terms proportional to the tree amplitude on the first line of eq. (8.9) get shifted 
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only in the denominator, namely only in S2...r in eq. (A9). This means that the overlap 
contribution from this term is the negative of the term itself and thus cancels it in the full 
amplitude, eq. (5.29). Furthermore, the same happens with the term involving h ^^. Thus, 
for these two contributions we have 

O'r = -(^ + 0Ar(l-,2-,3-,4+,...,n+), (8.12) 

ot-^- = I y: kt' — '"-\ ^^-^^^^^ — ^ • (8. 13) 

The d ^^-terra gives a contribution to the overlap, 

Qi-te™ = _1 V (3l^3....fe|l+) (3-|^3....^2....|l+)' (r(r + l)) ^ ^^_^^^ 

6 ^^^ (2+1 ^2...r k+) (2+1 ^2...r I (^ + 1) + ) f[ {k {k + 1)) S2....S3..,, 

k=3 

In discussing the ^f^^^'^-term, we first re-index eq. (8.3) by letting j ^ j — 3, so the sum 
runs from j = 4 to j = r. The ^f^^^'^-term has three different kinds of poles, since S2...j, S2...r 
and S2...(r+i) are shifted, 

;,W = - J-"^'' (8.15) 

(1 I -^2...(r+l) |2 ) 

corresponding to the overlap diagrams in fig. 17. At first glance it may seem that for j = r 
we obtain a double pole, since the denominator contains S2 r- However, from the numerator 
of the last factor in ^ , we obtain a factor S2...r, which cancels one of the factors in the 
denominator, and so we have again only single poles. Nevertheless, we treat the j = r term 
separately and obtain for this contribution, 

"-2 /o-l r/- r/ M-l-\ /o-l r/- r/- M-l-\2 
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fc=3 



(8.16) 



After some straightforward algebra, we find the remaining contributions, 

^^^ ijij + l)) 



o^' - ,EE 



.=5 ,=4 (2+1 ^2..., |J + ) (2+1 $2..., \{J + 1) + ) n (M^ + 1)) S3..., S2.. 
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- (3-| (^,^2...r - t^..Jr) |l + > (1^1 ^2...r^0-+l)...r^2..., |2'> \ • 

In term (c) we have relabeled r + 1 -^ r in the sum. 

Combining the various overlap contributions from eqs. (8.12)-(8.14) and (8.16)-(8.19) 
gives us the complete overlap contribution, 

0„ = O^"'' + O^*^"" + Of-*^"" + 0(f)'(^)'^=^ + 0(f) + Oi^) + Oi=) . (8.20) 



D. Assembling the Three-Negative Helicity AmpHtude 

The full amplitude A^^°(l^, 2~, 3^, 4+, . . . , n"*") is obtained by combining the pieces ac- 
cording to eq. (5.29). In this equation, the completed-cut terms and their large-z behavior, 
C„(0) and Inf [i_2) O^, respectively, are given in eq. (8.1) and eq. (8.10). The recursive con- 
tribution R^ is given by eq. (8.11), and the overlap contribution 0„ is given in eq. (8.20). 
Finally, the large-z contribution Inf[i 2} A^-i^ is given in eq. (7.4). 

E. Four-Negative Hehcity Amplitude 

The tree amplitudes for this configuration were first computed in ref. [68], along with all 
split helicity configurations, A„(l^, 2^, . . . , m~, {m + 1)+, . . . , n+), where legs of like helicity 
are nearest neighbors in the color ordering. At one loop, the coefficients of all logarithmic 
terms in the split helicity amplitudes were then computed in ref. [26]. Thus, we have the 
completed-cut terms of the amplitude with four negative helicities. It is straightforward to 
subtract the spurious large- 2; behavior of the completed-cut terms under a [1, 2) shift, and to 
extract the overlap terms. The terms with large-z behavior under the [1, 2) shift are given 
in eq. (7.14), leaving only the direct recursive contributions to be computed, according to 
eq. (5.29). 

The [1,2) shift recursion relation for four color-adjacent negative-helicity gluons may be 
read off from fig. 16, 

i?^(l-,2-,3-,4-,5+,...,n+) = 

n-l 

Y, AT^,^^{r, i^^:..., {^ + 1)+, • • • , n+)^i?.(2-, 3", 4', 5+, . . . , ^+, -KtJ 
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n-2 



n-2 



, . . . , 
X Af-(2-,3-,4+,5+, . . . ,^+, -it^:.,) , (8.21) 

where we omitted vanishing diagrams. We have solved this recursion relation through n = 
8. Although we will not present the analytic solution here, in section IX we present the 
numerical value of Agl^^'^{l^, 2^, 3^, 4^, 5+, 6+, 7+, 8+) at one phase-space point. 

F. Consistency of the Computed Amplitudes 

In order to confirm our procedure for constructing the amplitudes, we have performed a 
number of checks on the amplitudes. A simple check is on the reflection symmetry of the 
amplitudes under a reversal of legs, 

An;i{l-, 2", 3", 4+, . . . , n+) = (-1)"+M„;i(3-, 2", 1", n+, . . . , 4+) , (8.22) 

and 

An-iir, 2-, 3-, 4-, 5+, . . . , n+) = (-1)"+M„;i(4-, 3", 2", r, n+, . . . , 5+) . (8.23) 

Our choice of shifts makes this symmetry unobvious, so confirming the symmetry provides 
a rather non-trivial consistency check. We have confirmed the above symmetry numerically 
through n = 10 for the case with three negative helicities and at n = 8 for the case with four 
negative helicities. We have also confirmed numerically that the amplitudes have the correct 
factorization properties with real momenta through n = 8. We checked that all unphysical 
spurious poles are removable, for n = 6 (as mentioned in section IV) and for n = 7. These 
checks provide a strong confirmation that we found the correct analytic expressions for the 
amplitudes. Moreover, they also validate our recursive determination of the amplitudes' 
large-z behavior. 
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IX. NUMERICAL POINTS IN PHASE SPACE 

In order to aid the future implementation of these amphtudes in numerical codes, we 
present values of the amplitudes at one point in phase space. We present numerics for 
selected helicity amplitudes up to eight points, as a demonstration of the utility of our 
methods for high-multiplicity processes. At six points, where independent numerical values 
are available [3], we can use these as an additional check on the amplitudes. 

We quote the numerical results in an unrenormalized form. The renormalization amounts 
to carrying out the subtraction of eq. (2.14). 

Since we have exact analytical expressions for the amplitudes, it is a simple matter to 
evaluate the amplitudes to arbitrary precision; we give ten significant digits. High precision 
can be useful for studying the properties of the amplitudes near (removable) spurious singu- 
larities, in order to investigate numerical instabilities (roundoff error, etc.) which might be 
encountered under phase-space integration. Because our analytic expressions possess only 
a relatively mild set of spurious singularities, in comparison to more direct evaluations of 
Feynman diagrams, we do not anticipate any significant complications arising from roundoff 
error when constructing an NLO program. 

At six points, for ease of comparison, we choose the same numerical point as given in 
ref. [3]. 

— (—1, sin^, cos ^ sin 0, cos ^ cos 0) , 

— (—1, — sin^, — cos^sin0, — cos^cos0), 
I (1,1,0,0), 
-(1, COS/9, sin/3, 0), 

— (1, cos a cos (3, cos a sin /3, sin a) , 
6 

-ki - A;2 - /cs - ^4 - ^5 , (9.1) 

where 

9=-, 0=-, a = -, cosf3 = -—. (9.2) 

4 ' ^ 6 3 ' '^ 19 ^ ' 

With this choice, the energies of ki and k2 are negative, representing a physical scattering 

process at a collider. As in ref. [3] , we choose fi = n = 6 GeV, where /i is the scale originating 

from the dimensionally-regulated integrals. 
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ki 
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k2 
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ks 
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k. 
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k. 


= 


ke 


= 


1 = 


4 



At seven points, we choose the kinematic point, 



ki = — {—1 , sin 6 , cos 6 sin (f) , cos 6 cos (f)) , 

k2 = — (— 1, — sin^, — cos^sin0, — cos^cos< 

ks = |(1, 1,0,0), 

/c4 = — (1, cos/3, sin/3,0) , 
8 

/c5 = — (1, cosq;cos/3, cosQ;sin/3,sinQ;) , 

/ce = — (1, cos 7 cos /3, cos 7 sin /3, sin 7) , 



ky = —ki — k2 — k^ — k4 — k^ — ke , (9.3) 



where 



/I ^ J ^ ^ 27r ^37 ,^ ,, 

6 = —, (h = — , a = — , 7 = — , cos/3 = , (9.4) 

4 ' ^ 6 ' 3 ' ' 3 ' '^ 128 ' ^ ' 

and fi = 7 GeV. 

At eight points we choose as our reference kinematic point ^, 

ki = — (— 1, sin^, cos^sin0, cos^cos0) , 

/c2 = —(—1, — sin^, — cos^sin0, — cos 6* cos 0), 

k, = |(1, 1, 0, 0), 

A;4 = —(1, cos/3, sin/3, 0) , 
o 

/cs = — (1, cosa;cos/3, cosQ;sin/3, sina) , 

/cg = —(1, cos 7 cos /3, cos 7 sin /3, sin 7) , 
5 

kj = —(1, cos 5 cos /3, cos 5 sin /3, sin 5) , 
6 

kg = -ki - k2- ks- k4- k^- ke- kr, (9.5) 

where 

^ = ^, '^^e' ""3' ^^Y' — 3"' cos/3 = -— , (9.6) 

and we choose /i = 8 GeV. 



^ For this kinematic point, the quantities (7^|(1 + 2)|5^) happen to vanish; however, this vanishing causes 
no problems in evaluating the amplitudes considered here. 
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TABLE III: Numerical results for the nonvanishing A^ = 4 six- and seven-gluon split helicity 
amplitudes in the FDH scheme. The kinematic point is given in eqs. (9.1) and (9.3). The analytical 
expressions used for this table are from refs. [5, 6, 22]. 



Helicity 


lA 


eO 


— ++++ 


448.1350970 + i 288.8591589 


231.6837670 + i 1219.687214 


+++ 


- 12.83149626 + i 36.01649473 


- 82.74583978 + i 33.53625588 


—+++++ 


2923.502435 + i 683.4723607 


5112.775012 +i 6035.881921 


++++ 


- 45.77174817 + i 14.37948381 


- 118.6756853 - i 49.76271406 



TABLE IV: Numerical results for the nonvanishing Af = 1 chiral contributions to six-, seven- and 
eight-gluon split helicity amplitudes. The kinematic point is given in eqs. (9.1), (9.3) and (9.5). 
The analytical expressions were obtained from refs. [6, 24, 26, 49]. 



Helicity 


lA 


eO 


— ++++ 


- 28.11035867 + i 4.643367883 


- 108.9419206 + i 35.02980993 


+++ 


- 0.7843602985 - i 1.886592441 


0.1435789412 - i 5.332008939 


+++++ 


- 104.5611840 + i 45.34709475 


- 429.6932951 + i 209.0560823 


++++ 


1.069354012 -i 1.474267821 


1.619538641 + i 3.617458324 


—++++++ 


- 0.1599026654 - i 0.1215536994 


- 0.4460113565 - i 0.2626240110 


+++++ 


-0.01086792081 - i 0.002373299774 


- 0.003818661190 - i 0.02634017264 


++++ 


- 0.006820678758 - i 0.003628936145 


0.02086909214 - i 0.04355499656 



As explained in section II, it is convenient to decompose QCD amplitudes in terms 
of A/" = 4, A/" = 1, and A/" = amplitudes. QCD amplitudes are recovered from these 
components using eq. (2.13). 

We collect the numerical values of the split helicity amplitudes in Tables III-V. We have 
extracted an overall factor of icr from the numerical values presented in the tables. The 
complete set of analytic expressions for six- and seven-point A/" = 4 amplitudes was obtained 
in refs. [5, 6, 14, 22]. The A/" = 1 amplitudes appearing in Table IV were computed in 
refs. [6, 18, 24, 26, 49]. The finite A/" = amplitudes were obtained in refs. [66, 67]; a compact 
representation of the amplitude with a single negative-helicity was given more recently [31]. 
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TABLE V: Numerical results for the A^ = scalar contributions to six-, seven- and eight-gluon split 
helicity amplitudes. The analytical expressions were obtained from refs. [6, 26, 27, 31, 32, 66, 67] 
as well as sections IV and VIII of the present paper. 



Helicity 



lA 



++++++ 



0.1024706290 + i 0.5198025397 



-+++++ 



2.749806130 + i 1.750985849 



-++++ 



9.370119558 + i 1.547789294 



45.80779561 + i 13.03695870 



-+++ 



0.2614534328 - i 0.6288641470 



0.3883482043 - i 5.830791857 



0.1815778027 + i 1.941357266 



-++++++ 



22.52927821 + i 5.464377788 



34.85372799 + i 15.11569825 



176.2169235 + i 87.93931019 



-++++ 



0.3564513374 - i 0.4914226070 



0.7087164424 - i 11.32916632 



0.0009856214410 + i 0.002143695505 



0.001078316199 + i 0.03129931739 



0.05330088846 - i 0.04051789981 



0.05513350697 + i 0.1659518861 



-+++++ 



-0.003622640270 - i 0.0007910999246 



0.02719752089 - i 0.02586206549 



-++++ 



0.002273559586 - i 0.001209645382 



0.01154855076 - i 0.0008935357840 



The logarithmic parts of the divergent amphtudes were computed in refs. [6, 26], while the 
rational-function parts were determined in refs. [27, 32] and in sections IV and VIII of the 
present paper. Multiple analytical calculations exist for some of the split-helicity amplitudes 
appearing in the tables. In particular, the A/" = 4, A/" = 1, and logarithmic parts of the A/" = 
MHV amplitudes were also computed in refs. [11, 12] using MHV vertices [9]. 

Table III gives the numerical values for the 1/e and finite terms of the A/" = 4 super- Yang- 
Mills theory split-helicity six- and seven-point amplitudes. We do not include the coefficients 
of the leading 1/e^ singularities in the tables, as these are easily extracted from the values 
of tree amplitudes, for any helicity configuration. 



A/'=4 






^•^r .tree 



-Al 



(9.7) 



l/e2 e^ ■' 

The numerical values of the tree amplitudes may be read off from the values of the 1/e 
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singularities of either the A/" = 1 or the J\f = loop amplitudes, 



Ar=i 






cr 



cr 



-AT' , <r° , = ^^-"^^ , (9-8) 

l/e Ot 



lA e 

given in Tables IV and V. 

We quote numerical values of the Af = 1 and Af = amplitudes through eight points, but 
refrain from giving numerical values for eight-point A/" = 4 amplitudes. (The A/" = 4 ampli- 
tude with four negative helicities has not been evaluated explicitly, although it is straight- 
forward to obtain using the methods of ref. [16]. The general A/" = 4 n-point amplitudes 
with three negative-helicity legs may be found in ref. [37].) 

We have compared our numerical results for the six-point amplitudes to those found in 
ref. [3]. After accounting for differing overall phase conventions, they agree to within the 
number of digits quoted in ref. [3] . 

X. CONCLUSIONS AND OUTLOOK 

In this paper we have continued the development of the on-shell unitarity-bootstrap 
method [27] for computing complete amplitudes in non-supersymmetric gauge theories. It 
combines unitarity [5, 6] with on-shell recursion [28, 29]. It thereby systematizes an earlier 
unitarity-factorization bootstrap [20] used to compute the one-loop amplitudes required for 
Z —^ 4 jets and pp ^ W, Z + 2 jets. In the combined approach, cut-containing pieces are 
obtained via unitarity, while purely rational terms are obtained via on-shell recursion. The 
latter terms had previously been the most difficult part of a one-loop QCD calculation. The 
use of an on-shell recursion reduces them to tree-like calculations. 

On-shell recursion relations rely on continuations of the amplitudes to complex mo- 
menta [29]. As discussed in refs. [27, 30, 31], the complex factorization properties of loop 
amplitudes are much more subtle than those of tree-level ones. In addition to the presence of 
branch cuts, in loop amplitudes we encounter double and 'unreal' poles, present only when 
the momenta are complex. Unlike the situation for real momenta, there are as yet no gen- 
eral theorems describing how loop amplitudes factorize for complex momenta. For general 
helicity configurations, some continuations of the amplitudes may require such unknown 
factorizations. On the other hand, in general it is not possible to choose a continuation 
where all factorizations are known without spoiling the vanishing of the contour integral at 
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infinity whose consideration gives rise to the recursion relation. This unhappy situation can 
come about because the continued amphtude may not vanish as the continuation parameter 
becomes large. 

How, then should we avoid channels with unknown factorization, without spoiling the 
recursion? In this paper we provided a general strategy for avoiding these difficulties. We 
construct an auxiliary recursion relation to fill in information required to re-establish the 
contour integral and hence the primary recursion. The latter is designed to avoid channels 
with non-standard factorizations, even at the price of bad behavior as the continuation 
parameter becomes large. The auxiliary recursion supplies the terms we need to subtract 
from the contour integrand in order to obtain a well-defined integral and thereby a proper 
recursion relation. It may have channels with non-standard factorization, but is designed so 
that these non-standard channels will not affect the required subtraction terms. 

As an illustration of these ideas, we obtained complete n-gluon amplitudes with a 
scalar in the loop and three nearest-neighbor negative helicities in the color ordering. 
The logarithmic parts of these amplitudes were obtained in refs. [6, 26]; here we con- 
structed a recursive expression for all the rational terms that appear in these ampli- 
tudes. For the six-point scalar loop amplitude /lg^i"°(l~, 2~, 3~,4+, 5+, 6+) we found a 
compact representation for the amplitude. Using the supersymmetric decomposition [48], 
together with the previously computed A/" = 4 and A/" = 1 supersymmetric contribu- 
tions [6, 14, 22, 24, 37, 49], these provide complete QCD amplitudes. As an additional 
illustration, we provided a recursive expression for the eight-point amplitude with four- 
negative helicities, A^]"°(l^, 2", 3~,4~, 5+, 6+, 7+, 8+). We gave numerical values for the 
amplitudes (in the supersymmetric decomposition) at specific points in phase space, to use 
as a guide for constructing numerical programs in the future. In a companion paper we shall 
present all remaining n-point MHV amplitudes [47]. These calculations of high-multiplicity 
and indeed all-n amplitudes illustrate a crucial feature of our approach: the complexity of 
the calculation grows only modestly as the number of external particles increases. 

The continuations typically used to derive recursion relations are implemented by shifting 
a pair of momenta associated with external legs. We choose pairs of such shifts. Our strategy 
for choosing them is based on a several empirical properties of amplitudes: that under a 
[— , +) shift, gluon amplitudes are well behaved for large shift parameters, that the large- 
parameter behavior is always polynomial in the shift parameter, and that we can suppress 
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unwanted channels with non-standard complex factorization through appropriate choice of 
shifts. As yet we cannot offer formal proofs of these properties. In lieu of such proofs, we have 
checked the computed amplitudes by verifying the stringent real-momentum factorization 
properties any amplitude must satisfy. The correctness of the computed amplitudes may 
then be seen as further evidence for the correctness of the empirical properties. 

It would be useful to address some of these open formal issues, and to develop a first- 
principles understanding of the complex factorization properties and large-shift-parameter 
behavior of one-loop scattering amplitudes. Recent papers connecting tree-level on-shell 
recursion to the gauge theory Lagrangian [62] provide one possible avenue to deriving the 
analytic properties. Unitarity in D dimensions may also assist in this formal understand- 
ing [7], as it can be used to determine rational parts of amplitudes, even though explicit 
computations are often cumbersome. 

The methods we have presented in this paper are systematic and thus lend themselves to 
automation, most obviously in an analytic approach, but also plausibly in a semi-numerical 
one. The large number of required subprocesses in many applications of phenomenological 
interest makes such automation desirable. The techniques of this paper require as input the 
cut-containing parts of the amplitudes. Efficient, semi- automated procedures for evaluating 
unitarity cuts are therefore also necessary. As remarked in the introduction, there have 
been important developments in this direction in the past year or so, which are applicable 
to generic helicity configurations and external states in non-sup ersymmetric gauge theories, 
and which reduce the problem to one of residue extraction [18, 19]. In particular, this 
approach has been used to complete the computation of all cut-containing parts of six-gluon 
amplitudes. These methods should work well in tandem with the ones presented here. 

The empirical properties of the amplitudes are robust enough, and sufficiently broad, to 
allow the computation of arbitrary gluonic amplitudes. The same ideas should carry over to 
amplitudes with external quarks, vector bosons or Higgs particles. Beyond that, we foresee 
that with suitable extensions the method will also work for processes with massive particles 
propagating inside the loops. We are encouraged by the observation that the unitarity 
method applies to massive loops [7] and that at tree level, the on-shell recursive approach 
carries over to massive theories without difficulty [44, 45]. This extension of the on-shell 
bootstrap, to loop-level processes containing top quarks for example, is an open problem. 

Even before such extensions to amplitudes with massive internal particles become avail- 
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able, however, the on-shell bootstrap method we have described is ready to tackle many of 
the important multi-parton next-to-leading order computations of phenomenological interest 
at the Large Hadron Collider. 
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APPENDIX A: PREVIOUSLY COMPUTED AMPLITUDES AND VERTICES 



As discussed in the text, our bootstrap approach relies on using previously computed 
amplitudes. In this appendix we list the amplitudes that enter into our calculation. In 
the text, we also need the parity conjugates of the listed amplitudes, which are given by 
replacing {i j) *-^ [ji]. 



An{l^\2^\...,n^- 



Some useful tree amplitudes are 

Af"(r,2-,3+ 



Af%rXX 



Af<=(r,2-,3+,4+ 

Af^(r,2-,3+,4+,5+ 
Ar(l-,2-,3-,4+,5+ 



/lf^(r,2",3-,4+,5+,6^ 



[A„(l-'^^2-'^^...,n-'^^ 



(12r 



{i3)^\ji\ 



'(12) (23) (31)' 
[12] [2 3] [31]' 

'(12) (2 3) (3 4) (41) ' 

(12)^ 
'(12) (23) (34) (45) (51)' 

[12] [2 3] [3 4] [4 5] [51]' 

. (l'|(2 + 3)|4-)^ 

's234[23][34](56)(61)(5-|(3 + 4)|2-) 



+ i- 



(3-|(4 + 5)|6')' 



(Al) 



(A2) 
(A3) 
(A4) 
(A5) 
(A6) 



(A7) 



S345[21][16](3 4)(4 5)(5-|(3 + 4)|2-)' 

(The three-point amplitudes are non-vanishing for complex momenta, even though they 
vanish for real momenta.) The n-point MHV amplitudes [33, 69], and one sequence of 
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NMHV amplitudes [43], are, 

Ar(l^...,mr,...,m2-,...,n+) = ^^^^^^i^^^^ , (A8) 

UiHk + l)) 
fc=i 

Ar(l-,2-,3-,4+,...,n+) = (A9) 

(3-|^3....^2....|l+)'(r(r + l)) 



E 



r=4 (2+1 ^3..., |r+) (2+1 ^3..., |(r + 1) + ) S2...r S3...r R (^ (^ + 1)) 

fc=3 

Tree amplitudes with all like-helicity gluons or only one gluon of opposite helicity vanish. 

We will also need certain classes of one-loop helicity amplitudes, or rather, their rational 
parts Rn, defined in eq. (5.12). The following three- vertices with a scalar in the loop vanish, 

i?3(r,2+,-A1) = o, (AlO) 

related to the vanishing of the corresponding splitting amplitudes for real momenta [5]. 
In principle, there are also vertices with like-helicity external legs, -R3(l^, 2~, — P^2) ^^^ 
i?3(l+, 2+, — P^2)- However, the associated complex factorization properties are not fully 
understood as yet, so we avoid these channels, or suppress them in a large-z limit. 

For multi-particle factorizations, factorization functions can appear, as discussed in sec- 
tion V. For the scalar loop, the one-loop factorization function is, 

^(-'-(^5'"(^)-D' <-^' 

The two-point rational "vertex" associated with this is a constant. 

As discussed in section II, it is convenient to apply a supersymmetric decomposition of 
QCD loop amplitudes into A/" = 4, A/" = 1, and A/" = (scalar loop) pieces [48, 65]. Just as 
the tree-level amplitudes with all gluons of like-helicity or all but one gluons of like-helicity 
vanish, so do the corresponding supersymmetric A/" = 4 and A/" = 1 amplitudes; the scalar 
A/" = amplitudes are nonvanishing, but purely rational and finite. For the finite four-point 
helicity amplitudes we have [57, 65], 



<r°(l-,2+,3+,4- 



" 3 (12) (34)' 
. cr (2 4) [2 4]^ 
^ 3 [12] (2 3) (3 4) [41 ' 
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(A13) 
(A14) 



The corresponding five-point amplitudes are [48], 



4jV=0fT+ 0+ Q + 



A^='il\2\3\A\5' 



<r°(l-,2+,3+,4+,5^ 



.cr 



[1 2] [2 3] 



[4 5] [51] 



[2 5] [3 4] 



_(34)(45)(51) (12) (23) (34) (1 2) (34) (5 1) J ' 



(A15) 



cr 1 



3 (3 4)^ 



[2 5] = 



(14)^ [4 5] (3 5) (13)^ [3 2] (4 2) 



[12] [51] (12) (2 3) (4 5)' (1 5) (5 4) (3 2)'. 

(A16) 



Expressions for all remaining one-loop finite amplitudes may be found in refs. [30, 66, 67]. 
We make use of the following functions [48] to express the other, cut-containing ampli- 
tudes at loop level, 

,2 



Ko(s) 
Lo(r) 
Li(r) 
Uir) 



1 / 




e(l-2e)V 


ln(r) 




1-r' 




Lo(r) + 1 




1-r ' 




ln(r) - (r - 


l/r)/2 



\ + ^^(r^) + 2 + 0{e), 



(A17) 



(l-r)3 

in order to eliminate spurious singularities for r ^ 1 which are present in the pure cut terms. 
The four-point amplitude ^4.1(1^, 2^, 3+, 4+) is given by. 



Aff4(l-,2-,3+,4+) = crA*'^<='(l',2-,3+,4^ 



/i^ 



-S12 



+ 



/i^ 



V-S23 



"S23 



(A18) 

(A19) 

, ■ Q,, (A20) 

-S23/ 9/ 

where the regularization-scheme parameter 5r was discussed in section II. The recursive 
vertex we need for the A/" = scalar loop contribution is therefore given by 



<fi(l-,2-,3+,4+) = crA*'--(l-,2-,3+,4+)f- + lnf-^)+2), 

^e ^-S23^ ^ 

<r°(l-,2-,3+,4+) = crA*-^(l-,2-,3+,4+)(l + lln(- 



i?4(l-,2-,3+,4+)=(- + |)A*-<=(l-,2-,3+,4+). 
The remaining five-point amplitude with adjacent like helicities that we use is. 



(A21) 



A^f^(l-,2-,3-,4+,5+) = crA*^<=^(l-,2-,3-,4+,5 



+\ 
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X 



E 



li' 



<r'(l-,2-,3-,4+,5+) = 









In 



-■Si+2j-2 
-Sj-2,i-l 



5 9 5r 



-TT 



/i^ 



2e l V-S51 



+ 



/i- 



-S34 



, (A22) 



(A23) 



—I 



ct 



[45]^ [51] (12) [24] + [5 2] (23) [34] Lo'^ 



S34 



2 [1 2] [2 3] [3 4] [5 1] 

<r°(r,2-,3-,4+,5+) = 
\N=l 



S34 



-^^1=1(1-, 2-, 3-, 4+, 5+) + -crA*-^(l-, 2-, 3", 4+, 5 



+ i 



cr 



9 

(12) [2 4] [5 2] (2 3) f [5 1] (12) [2 4] + [5 2] (2 3) [3 4]) L 



(A24) 



-^51 
-S34 



[1 2] [2 3] 



cr-Rs, 



'34 



where 
i?5 = i 



[13](13)= 



(13)' [4 5] 



1(1 2) (2 3) [2 4] [2 5] [4 5] 



3 [12] [2 3] (3 4) (4 5) (51) 3 [1 2] [2 3] (3 4) (5 1) 6 [1 2][2 3]s34S5i 
The vertex we need is given by setting the logarithms to zero in the amphtude, 



i?5(l-,2-,3-,4+,5+) = -<r°(l",2-,3-,4+,5^ 

Cr 



ln=0 



(A25) 



(A26) 



The A/" = 4 supersymmetric amphtude needed for constructing the complete six-point 
QCD amplitude with three color adjacent negative helicities is given by [6], 



lA/'=4/i- o- Q- 



(1) 



.(2) 



A^=\r, 2-, 3-, 4+, 5+, 6+) = cr [B^ W'^'> + B^ W'^'> + B^ W, 



.(3)1 



where 



Bi = i 
B2 = i 



^■5123 J 



[1 2] [23] (45) (56) (4^1 (2 + 3) |r) (6-| (1 + 2) [3') 

-\3 



(A27) 



(A28) 



(11(2 + 3)14-) 



5234 [2 3] [3 4] (5 6) (61) (5- 1(3 + 4) |2-) 
(2 3)^5 6]^ 



+ t 



S234 (34) [6 1] (4-| (2 + 3) |1-) (2-| (3 + 4) |5-) 



(A29) 



5, 



(3-|(l + 2)|6-)' 



+t 



S345[61][12](3 4)(4 5)(5-|(6 + l)|2-) 
(12)' [4 5]' 



(A30) 



S345 (6 1) [3 4] (6-| (1 + 2) |3-) (2-| (6 + 1) |5-) ' 
We have taken the parity conjugate and somewhat simplified the results compared to the 
expression of ref. [6], and used the four- dimensional helicity scheme where 6^ = 0. The 
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combination of integral functions appearing in the amplitude is, 



^f —^e(^)' -M^i^y 



i=i 



2e ^-^ V — S^('^_L1^ / V —Sid-LT) 1 V —s 



i(i+i)/ V -^i(i+\) J V--^(j+i)(i+2) 



In ( ^fM!±^ iJ^fM!±^ + ln( ~''-^^^ In ^ "'^■■•(^+^) 



^S(i+3)(i+4)/ \— S(j+4)(i+5)/ \ — S(j+2)(i+3) / \^S(J+5)i 



2 \— S(j+3)(i+4)/ y— S(i+4)(i+5)/ 2 V^-^JCi+l)/ \~S(J4.2)(i+3) 

+ i In fz£(!±^)(!±2) ") In f Z£(!±l)(!±^') + t. . (A31) 

2 \-S(j+3)(i+4)/ V -S(i+5)j y 3 
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